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DEFINITION AND APPLICATION OF LONGITUDINAL STABILITY 
DERIVATIVES FOR ELASTIC AIRPLANES 

By William B. Kemp, Jr. 

Langley Research Center 

SUMMARY 

A review of past practices in the analysis of longitudinal stability of elastic air- 
planes has revealed some inconsistencies arising from an inappropriate constraint on 
airplane speed. These inconsistencies lead to ambiguous definitions of stability deriva- 
tives and inaccurate prediction of some stability characteristics and, possibly, even per- 
formance characteristics. 

A single, consistent set of longitudinal stability derivatives for elastic airplanes is 
defined in this paper. Unlike past practice, the aeroelastic contributions of dynamic- 
pressure perturbations are included, and the aeroelastic contributions of normal acceler- 
ation appear primarily in the derivatives with respect to pitching velocity and angle-of- 
attack rate. 

Approximate expressions for static stability and control parameters in terms of 
these derivatives are shown to correlate well with results of more complete dynamic 
solutions. Some results of an illustrative analysis cast doubt on the general existence of 
a neutral point or maneuver point and indicate the possibility of a strongly divergent phu- 
goid oscillation at relatively high frequency, particvilarly when the variation of atmo- 
speric properties with altitude are considered. One commonly used alternative formula- 
tion of the stability derivatives is shown to yield significantly inaccurate predictions of 
not only the speed stability and phugoid characteristics, but also the short -period fre- 
quency and damping. 


INTRODUCTION 

The overall flight performance of many aircraft designs of current interest is 
strongly influenced by the effects of elastic deformation of the aircraft structure imder 
the action of the flight loads. These aeroelastic effects are particularly important for 
aircraft having large size, low structural design load factor, low structural weight frac- 
tion, and high operating dynamic pressure. Aeroelastic effects are manifested not only 
in the area of structural dynamics (flutter, landing loads, etc.) but also in the quasi-steady 
areas of stability and control and trimmed flight performance. 



Two different approaches are widely used for analysis of quasi-steady aeroelastic 
effects: the modal technique and the influence -coefficient technique. As pointed out in 
reference 1, both techniques can utilize the same description of aerodynamic, elastic, 
and mass properties but differ in the form of the perturbation variables whose effects 
are superimposed to describe the structural deformation. In the modal technique, a rela- 
tively complete dynamic formulation involving the normal modes of free vibration is 
reduced to a quasi-steady formulation by assuming that the velocities and accelerations 
in these vibration modes are negligible. The influence -coefficient technique applies 
steady-state solutions for structural deformation to the relatively low frequency rigid- 
body motions. Both techniques are developed in some depth in reference 2, which points 
out that from either technique, the quasi-steady (termed "equivalent elastic" in ref. 2) 
formulation embodies the implicit assumption that the structural deformation is in phase 
with the applied loads and therefore the dynamic analysis requires no more degrees of 
freedom than are needed for analysis of a rigid airplane. The quasi-steady effects of 
aeroelasticity are then represented completely by changes in the aerodynamic stability 
derivatives and by new derivatives which are either nonexistent or negligible for a rigid 
airplane. 

A review of the procedures outlined in references 1, 2, and 3 for calculating the 
longitudinal stability and control characteristics of elastic airplanes reveals some dis- 
turbing inconsistencies. For example, in reference 1 the forward-speed degree of free- 
dom is omitted from the basic dynamic formulation, yet a static stability parameter for 
Ig level flight is developed. This omission of the speed degree of freedom is believed 
to be typical of many applications of the modal technique. In reference 2, aerodynamic 
derivatives with respect to dynamic pressure are included in the expression for longitu- 
dinal static speed stability but are omitted in the dynamic stability analysis. Moreover, 
all the derivatives recommended for use in the dynamic analysis are derived imder a 
constraint which is valid only at constant dynamic pressure. Reference 3 makes use of 
the same constraint but suggests the use of two neutral -point locations, one applicable to 
constant -speed maneuvers and the other applicable to constant -load -factor maneuvers. 
These inconsistencies can lead to errors not only in the predicted stability and control 
characteristics, but also in the predicted trimmed drag coefficient through use of inap- 
propriate aerodynamic derivatives related to lift -curve slope and static stability. 

When it is considered that meaningful indices of static stability dvu'ing perturbations 
in either speed or load factor can be obtained from control response transfer fimctions 
resulting from a single solution of the dynamic longitudinal equations of motion, it is 
apparent that a single set of aerodynamic stability derivatives must exist for a given flight 
condition. These stability derivatives are then applicable, within the limitations of linear 
dynamic analysis and the quasi -steady aeroelastic assumptions, to the prediction of all the 
pertinent static and dynamic stability and control characteristics as well as aerodynamic 


performance characteristics. In this paper, an attempt to define such a set of derivatives 
is described and the significance of certain of the derivatives to the longitudinal stability 
and control problem is examined by means of illustrative calculations. 

SYMBOLS 


a 

[A] 


a 


[B] 




b 


Ca 


Cd 

Cm 

Cn 

Cx 


Cz 


c 

F 


aerodynamic operator 
aerodynamic influence -coefficient matrix 
velocity of sound, meters/second 
aeroelastic correction matrix 

matrix used in calculating dynamic -pres sure contributions to stability 
derivatives 

matrix used in calculating Mach number contributions to stability derivatives 


reference span, meters 
axial -force coefficient, 
drag coefficient, 


•X,A 


lift coefficient. 


Lift 

qS 


pitching -moment coefficient, 
normal -force coefficient, -- 
local X-force coefficient, 
local Z -force coefficient, 
reference chord, meters 
force, newtons 


My, A 
qSc 
'’z,A 


qS 

^X,A 

qS^ 

^Z,A 

q^ 
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Fx»Fy,Fz 


force components along reference axes, newtons 


g gravity acceleration, meters/second^ 

gjj gravity acceleration at sea level, meters/second 

H angular momentum, kilogram -meter s^/second 

h altitude, meters 

T inertia dyadic 

Ix,Iy>Iz moments of inertia about reference axes, kilogram-meters^ 

Ixz product of inertia about X- and Z-axes, kilogram -meter s^ 

T,T,k unit vectors along reference axes 

My + My 2 

M Mach number ; or 2—^^ 1/second 

Mx?My,Mz moment components about reference axes, newton-meters 
m mass, kilograms; or any motion variable 

n . component of combined kinetic and gravity acceleration in negative 

Z -direction, normalized by gg 

p,q,r angular velocities about reference axes, radians/second 
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XI O 


J 


structural operator 


[S] 


t 

u,v,w 

u 


structural slope influence -coefficient matrix 
Laplace variable 

displacement of propulsive thrust control 
time, seconds 

velocity components along reference axes, meters/second 

V - Vi 


normalized velocity perturbation. 


velocity, meter s/second 


Vl 


w 


airplane weight, newtons 


X,Y,Z 


orthogonal reference axes 


Fy A + Fy rri 

X = 1/second 


mV 

M 

x,y,z coordinates relative to reference axes, meters 

distance behind leading edge of reference chord, fraction of reference chord 

F7 A + F7 m 

Z = l/second 

mVj^ ’ 


a 

/3 

6 


^6 


angle of attack, radians 
angle of sideslip, radians 

angular displacement of aerodynamic control surface, radians 
local slope of mean camber surface, - ^ 

weighted value of mean -camber -surface slope used to represent unit 6 
damping ratio of phugoid mode 
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^sp 

M 

P 

a 

u> 

Wp 

^sp 

Subscripts: 

A 

a 

c.g. 

des 

F 

f 

I 

ij,k 

jig 

I 


damping ratio of short -period mode 

Euler attitude angles in pitch, roll, and yaw, respectively, taken in order 
of \{/, 9 , 4 > from earth reference to airplane reference, radians 

air viscosity, newton-seconds/meter^ 

air density, kilograms/meter^ 

induced downwash angle used as boundary condition, radians 
angular velocity, radians/second; or circular frequency, radians/second 
undamped natural frequency of phugoid mode, radians/second 
undamped natural frequency of short -period mode, radians/second 

from aerodynamic sources 
value at aerodynamic control point 
value at center of gravity 
value at airplane design condition 
increment due to flexibility 
value at force application point 
from inertial sources 
index indicating a member of a set 
value corresponding to jig shape 
local value 



p pertvirbation from reference flight condition 

S steady state 

s referred to stability axes 

T from propulsive sources 

U unsteady 

1 value at reference flight condition 

Mathematical notation: 

derivative with respect to time 
second derivative with respect to time 
vector quantity 
column matrix 
row matrix 
matrix 

diagonal matrix 

Stability derivative notation: 

Any of the following variables: 

# 

6,Q!,Qf,0,q,q,u,u,T,h 

when used as a subscript on any of the following quantities: 

X, Z, M, Cj^ Cfij, Cjn 

denotes the stability derivative of the indicated quantity with respect to the indicated 
variable. 


(*) 

n 

D 

o 

if 

[] 

[] 



DEVELOPMENT OF FORMULATION 


Fimdamental Assumptions 

As indicated in the introduction, it is desired to define a self-consistent set of aero- 
dynamic stability derivatives applicable to the calculation of the longitudinal static and 
dynamic stability and control and trimmed flight characteristics of an elastic airplane. 
These derivatives will be defined in the context of their role as constant coefficients in 
the dynamic equations of airplane motion. Since the aerodynamic processes involved are 
not necessarily linear with the motion variables, the concept of local linearization will be 
used so that the aerodynamic derivatives can be considered constant for small motion 
perturbations from a reference flight condition. The linearized small -perturbation equa- 
tions of motion then constitute an appropriate form for use of the derivatives. 

The basic assumption of quasi-steady aeroelasticity will be made, namely, the 
structural deflection is assumed to be in phase with the applied loads. This assumption 
is equivalent to ignoring the velocities and accelerations in the elastic degrees of free- 
dom. One result of this assumption, stated in reference 2, is that the set of equations of 
motion is reduced to include only those governing motion of the airplane center of gravity 
(abbreviated c.g.). The quasi -steady aeroelastic effects are then embodied completely 
in the aerodynamic characteristics, and the equations of motion are otherwise identical 
to those for a rigid airplane. The related assumption that the aerodynamic load change 
arising from structural deformation is in phase with the structural deflection is also 
made. In other words, the unsteady aerodynamic effects related to structural deforma- 
tion are neglected. However, this does not preclude the existence of aerodynamic deriva- 
tives with respect to the rate of change of angle of attack. 

Although use of the quasi-steady assumptions clearly introduces no error in the 
prediction of steady-state flight characteristics, some inaccuracy can occur in the char- 
acteristic stability (rigid body) modes because dynamic coupling between these modes and 
the elastic modes of structural vibration is neglected. This coupling can become increas- 
ingly significant as the frequency of any structural mode approaches that of any rigid-body 
mode. In spite of the possible inaccuracy, the quasi-steady assumptions are frequently 
made, particularly in preliminary aircraft design studies, because of the simplified 
dynamic formulation that results from their use. 

To describe the elastic properties of the airplane structure, a structural reference 
plane is assumed to exist and to be oriented in such a way that the in-plane flexibility is 
negligible relative to that normal to the reference plane. Furthermore, the longitudinal 
aerodynamic characteristics are assumed to be unaffected by elastic distortion of either 
vertical fins or fuselage cross sections. All remaining structural deflections are con- 
strained to the direction normal to a single plane chosen to lie as nearly as possible 
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parallel to the mean camber surfaces of the wing, body, and horizontal tail. It is appro- 
priate, then, to choose an orthogonal system of reference coordinate axes, as shown in fig- 
ure 1, so that the Z-axis is parallel to the direction of allowed structural deflections and 
the X-axis lies in the plane of symmetry. For small perturbations from the deflections 
existing in a reference flight condition, the center -of -gravity location, the pitching moment 
of inertia, and the aerodynamic and structural influence coefficients can be considered to 
be constant. In practice, these quantities, determined for a given reference condition, 
are frequently used over a set of flight conditions to be analyzed. 

Formulation of Stability Derivatives 

In airplane longitudinal stability analysis, one is concerned with the motion of the 
airplane as a whole in three degrees of freedom, including two components of translational 
freedom and one of pitch rotation. An additional degree of freedom exists for each inde- 
pendent control system. In the present development, two such control systems will be 
assumed, one aerodynamic and one propulsive, having control displacements denoted by 
6 and T, respectively. The conditions for dynamic equilibrium in the three airplane 
degrees of freedom are 


^X,A + %,T + ^X,l = 0 

(la) 

Fz,A + ^Z,T + Fz;j=0 

(lb) 

My,A + My,T + My, I = 0 

(Ic) 


where the nine force and moment components may vary with time. Conditions for dynamic 
equilibrium in the two control degrees of freedom are assumed to exist and be satisfied 
although they will not be specified. 

Equations of motion are formed from the equilibrimn equations when the forces and 
moments are expressed as ftmctions of a set of motion variables. The fundamental set of 
motion variables consists of the displacement and its successive time derivatives in each 
degree of freedom. Alternative sets of motion variables can be derived from the funda- 
mental set by transformation. The total number of motion variables is somewhat arbi- 
trary and depends on the maximum and minimum order of differentiation of each displace- 
ment needed to express the forces and moments with svtfficient accuracy for the problem 
at hand. Second derivatives of the airplane displacements are needed to define the inertial 
contributions and are believed to be of sufficiently high order to describe the aerodynamic 
and propulsive contributions at frequencies compatible with the quasi-steady assumption 
made in the present problem. At the low -order end, one or more of the displacements 
themselves may be omitted from the set of motion variables. Under the frequently made 
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assiunption of a uniform atmosphere, the forces and moments do not depend bn the loca- 
tion of the airplane, and even if atmospheric variations with altitude are considered, the 
horizontal displacement is of no significance. 

Alternative sets of motion variables can be derived by transformation from the 
fimdamental set. Common practice in the analysis of longitudinal aircraft motions has 
resulted in the use of variables such as (V,a, 6) as the set of lowest order motion vari- 
ables in the three airplane degrees of freedom. These variables must be supplemented 
by the control displacements 6 and T for the five -degree -of -freedom problem at 
hand. These five variables and their time derivatives constitute an appropriate set of 
motion variables only if a locally uniform atmosphere is assumed. If the variation of 
atmospheric properties with altitude is to participate in the dynamic problem, a more 
appropriate set will be (V,h,o!, 0,6,T), but because these variables are not mutually inde- 
pendent, the equations of motion must include a kinematic equation relating h to V, a, 
and 0. This scheme is retained in the present development. 

Common practice in dynamic stability analysis has led to the use of stability deriv- 
atives expressed in either coefficient form or a dimensional form having units of negative 
powers of time. The coefficient form is usually applied only to aerodynamic forces and 
moments, whereas the dimensional form includes both aerodynamic and thrust contribu- 
tions. To define these derivatives, equations (1) may be restated as 


^ ° + ^x,t) = 

^ + ^z,t) = - 

1 / \ I 


(2a) 

(2b) 

(2c) 


Expressing X, Z, and M as fxmctions of the motion variables and applying the prin- 
ciple of linear superposition of small perturbations give 


X^mj^ . . . mj . . . mjj = X^mj j . . . mj j . 
Z(mi . . . mj . . . mj) = Z^mi 1 . . . mj^j . 
M^mj i. . . mj . . . mjj= . . . mj j . 



(3a) 

(3b) 

(3c) 


10 



I 


where i is the value of the jth motion variable in the reference flight condition and 
mpj is the pertmrbation of the jth motion variable from its reference value. The partial 
derivatives appearing in equations (3) are defined as the complete set of stability deriva- 
tives in dimensional form. The stability derivatives in coefficient form are similarly 
defined as the partial derivatives of the coefficients Ca, Cfj, and Cm with respect to 
motion perturbations where 

(4a) 


Cn = - 


^Z,A 

qs 


, _ ®^y,a 


(4b) 

(4c) 


In determining the values of the stability derivatives by theoretical, experimental, or 
combined means, the set of motion variables is not necessarily the most appropriate set 
for isolating and understanding all the contributions to the stability derivatives. A more 
appropriate set might be those variables describing the physical processes involved in 
the production of aerodynamic and propulsive loads. In general, any change in the shape 
of an airplane surface will produce a charge in the aerodynamic load distribution. There- 
fore, any process capable of deforming the structure of an elastic airplane will contribute 
to one or more stability derivatives. Conversely, for an elastic airplane, a change in 
aerodynamic or thrust loading will, in general, change the structural deformation. Con- 
sequently, all stability derivatives reflect contributions from one or more processes 
capable of deforming the structure. Therefore, a set of physical variables capable of 
defining the structural deformation is a suitable set for determining the stability deriva- 
tives of an elastic airplane. 

If it is assumed that such a set of physical variables P|^ exists and that it can be 
functionally related to the set of motion variables, then a stability derivative, dX/amj 
for example, can be expressed as 



The number of physical variables K need not equal the number of motion variables J 
and the physical variables need not be mutually independent imder the constraints of the 
airplane dynamic problem. 

Under the concept described, the distribution of structural deflection can be consid- 
ered as the linear superposition of a number of deflection distributions, each identified 
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with a particiilar physical variable. For small perturbations from a given reference 
flight condition, the shape of each distribution is invariant and its amplitude is propor- 
tional to the perturbation of its physical variable. 

An appropriate set of physical variables may now be determined. Consider an air- 
plane configuration represented by a nearly planar mean camber surface divided into a 
large number of elemental panels. In symmetric flight, the configuration shape is defined 
for aerodynamic purposes by the paneling geometry and the set of chordwise slopes 
{e} = {-dz/ox} at all panels. The configuration is flying in an aerodynamic environ- 
ment defined globally by Mach number, Reynolds number, and dynamic pressure, and 
locally by the set of local angles of attack and angle -of -attack rates act- 

ing on each panel and arising from aircraft motion and induced effects of the propulsion 
system. The sets of X- and Z -components of the aerodynamic force coefficients acting 
on the panels can be expressed by 


{Cx,a} = ^X,S + ^X,U {^l} 

(6a) 

{Cz,a} = ^z,s («z + + ^Z,U 

(6b) 


where Cl represents an aerodynamic operator dependent on paneling geometry, Mach 
number, and Reynolds number. The subscripts S and U refer, respectively, to 
steady-state and first-order unsteady aerodynamic processes. The corresponding sets 
of force components are 


{Fx,a} = q {Cx,aSz} 

(7a) 

{Fz,a} = q ^Z,A^l} 

(7b) 


The distributions of local angle of attack, angle -of -attack rate, and surface slope 
are represented by the superposition of specifically identifiable contributions as follows: 

“ (l) - ^{xa> + (8a) 

{«Z> =«(l> -|{xa> (8b) 

("> = {^jig> + « ("6> + {^F> (8) 

The contributions to o!^ and arise from the time -dependent airplane angle of attack 

and pitch rate and the thrust -induced flow perturbations. The set {x^^ consists of the 
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x-coordinates of points (one per panel) at which the aerodynamic boundary condition of 
no flow through the surface is satisfied. The contributions to local surface slope, defined 
at the same points as Xa, arise from the jig shape (mean camber surface of the com- 
pletely unloaded structure), control -surface deflection, and structimal flexibility. 

The flexible contribution to surface slopes can be expressed as the influence of a 
structural operator jS operating on the Z -component of the forces acting on each panel 
arising from aerodynamic, inertial, and propulsive sources as follows: 

Z,A + Fz,I + F Z,t} (10) 

The aerodynamic -force set {Fz^a} been defined by equation (7b). The inertial - 
force set {F 2 results from the action on the mass distribution of the set of 
Z -components of linear acceleration arising from gravity and airplane motion. By using 
a lumped parameter concept, let mj be the mass assigned to the ith panel, assumed to 
be located at the point (xf^i,yf^i,0^. The effect of this mass can be replaced by the 
inertial -force vector 

•..■■",( 5 -$) 

where Vj^ is the velocity of m^ in inertial space. Under the quasi-steady assumption, 
the velocity and acceleration of mj[ relative to the airplane as a whole are neglected, 
and the Z -component of inertial force can be written (ref. 4) as 

Fz,I,i = cos e cos 0 - w - pv + qu + Xf^^(q - pr) - yf^|(p + qr)J 

In the following section, equations of motion applicable to analysis of longitudinal motions 
in a banked turn are developed. For this case, the bank angle cp and the rolling and 
yawing velocities p and r are nonzero but are considered to be constants. In equa- 
tion (30) the first four terms in the above equation are related to n, the normal acceler- 
ation sensed at the airplane center of gravity. The last term gives rise to an antisym- 
metric distribution of inertial force whose effect on the longitudinal stability derivatives 
should be negligible and will be neglected. With these considerations, the inertial -force 
set can be written as 

(Fz,]} = ngo [m] (l) + (q - pr) [n^ (xf) (11) 

where is the diagonal matrix of masses assigned to the individual panels and ^xj) 

is the set of x-coordinates of points at which loads are applied to the panels and is not 
necessarily identical with . The relationship of Xa and xf to panel geometry is. 
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of course, involved in the processes represented by both the aerodynamic and the struc- 
tural operators, Cl and jS . 

The propulsive -force set {Fz,t)’ Z -components of panel forces 

statically equivalent to the direct thrust reactions at the engine mounts. 

The quantities X, Z, and M and the coefficients C^, Cj^, and Cjjj can be 
obtained by summation over the set of panels lying on one side of the plane of symmetry. 
Thus, 


^ = 5^ (l)^(^X,A + %,T> 

(12a) 

0’z,A + ^Z.t) 

(12b) 


(12c) 

Ca = -i{Si)’'{Cx,A> 

(13a) 

Cn= -|<Si>’'{Cz,A} 

(13b) 

Cm = -^{%Jtl}^{Cz,A) 

(13c) 


It is now possible to identify the set of physical variables with which structural 
deflection modes can be associated and aerodynamic and thrust loads can be determined. 
Table I lists these variables and the symbols used for each value in a steady reference 

TABLE I. - THE PHYSICAL VARIABLES 


Physical variable 

Value in reference 
condition 

Perturbation from 
reference condition 

Equation 

Control deflection 

5l 


(9) 

Angle of attack 

“1 


(8a) 

Angle -of -attack rate 

0 

a 

(8b) 

Pitch rate 

‘ll 

% 

(8a) 

Pitch acceleration 

0 

• 

q 

(8b), (11) 

Normal acceleration 

«1 

np 

(11) 

Mach number 


Mp 

(6a), (6b) 

Reynolds number 

R1 

Rp 

(6a), (6b) 

Dynamic pressure 

^1 


(7a), (7b) 

Thrust command 

Tl 

Tp 

(8a), (10) 
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flight condition and for the perturbation of each from the reference condition. The equa- 
tions through which their effects enter the aeroelastic problem are also identified. 

In the foregoing discussion, the propulsive process was not described in detail. To 
define the form of stability derivatives, the assumption is made that the propulsion sys- 
tem outputs, {Fx,T^> ^Z,t} described in terms of the physical 

variables listed in table I. 

In concept, the partial derivatives of X, Z, and M or of Ca, Cn, and Cm 
with respect to each physical perturbation variable can be obtained through application 
of equations (6) to (13). The desired stability derivatives are then defined by relations 
such as equation (5) once the partial derivatives apj^y^amj have been determined. 

Because a mutually dependent set of motion variables has been adopted, the partial deriv- 
atives will be evaluated from expressions for each physical variable in terms of a mutu- 
ally independent set of motion variables obtained by omitting either h or a from the 
set (V,h,Q', 0,6, T). Thus, the physical variables 5, a, a, and T appear without change 
in the set of motion variables. The variables q and q will also appear unchanged in 
the motion variables, replacing 9 and 0 since, for longitudinal motions, 0 is a lin- 
ear fimction of q, as shown in the following section. The normal acceleration n is 
expressed as a function of V, a, and 0 and their time derivatives in equation (34) of 
the following section. The remaining physical variables are expressed as functions of 
V and h as follows: 


M = -^ 
a(h) 

(14) 

R _ p(h)Vc 
p(h) 

(15) 

^ _ P(h)v2 
q - 2 

(16) 

In applying the principle of small perturbations, 
dimensional perturbation in velocity magnitude as 

it is convenient to define a non- 

. V - Vi 

u =• 

Vi 

(17) 


where the perturbation velocity component uV^ is directed along the unperturbed veloc- 
ity vector and is orthogonal to the component oipVi, as shown in figure 1. 

After expressing each motion variable as the sum of a reference value and a per- 
turbation value, the partial derivatives /apj^ /am^V can be derived and are given in 
table n. ^ 
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TABLE n.- THE PARTIAL DERIVATIVES 



2qi U U U U ^qj^cosQ!j u u 

^ »o 

Vi 

0 0 0 0 0 - ^ sin Q! 1 0 0 

&o 

Vi 

0 1 0 0 0 -^qjSino!j 0 0 


0 0 1 0 0 - cos a I 0 0 


0 0 0 0 0 -^ sin 0i cos 4> 0 0 

°o ^ 

Vi 

0 0 0 1 0 ^ cos a I 0 0 


0 0 0 0 1 


^1 ^ 8a \ ^l/QP \ 0 0 0 0 

' ^lV9h/j Pi\9h/j " Pi\9h/j Pi\9h/j 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 


The desired stability derivatives can now be defined by using equation (5) and 
table n. In dimensional form, 


-#-sin 01 cos 





I 



where the subscript 1 denotes a quantity evaluated at the reference flight condition or, 
more specifically, at the reference condition values of the physical variables as noted in 
table I. The Z and M stability derivatives are defined by the expressions obtained 
by replacing X in equations (18) by Z or M. The defining expressions for the coef- 
ficient stability derivatives are equivalent to those for the dimensional derivatives but 
may include constant factors as appropriate to render each motion variable dimension- 
less; thus 


Ca. = 

2Vi/8Ca\ 


fdC^ 

c ^84 

Sq^ 

an y 

o 

11 

2Vi/aCA\ . 

2Vi^ 

onci *• 

c [dqj /- 

goC 

cos 


(19) 


The motion variables 6, a, d, u, and possibly T are already dimensionless. The 
stability derivatives of Ca, Cn, and Cm with respect to q, u, and h are considered 
herein to retain units reciprocal to those of the respective motion variables. 

In evaluating the stability derivatives from wind-tunnel experiments and/or theory, 
it is anticipated that the coefficient form will be evaluated first and then converted to the 
dimensional form for use in stability analysis. Expressions for the dimensional deriva- 
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tives in terms of the coefficient derivatives are needed therefore. Combining equa- 


tions (2) and (4) gives 



(20a) 

mVi mVj 

^ _ qs% Fz,t 

^ ~ mVi ' mVi 

(20b) 

qScC^ My^T 
ly ^ ly 

(20c) 


The partial derivatives with respect to the physical variables are 

M.= 9q qs i 

apjj mVj 9pjj mVj 9pjj mVj 9pj^ 

9 Z _ 9q qS 1 9Fz,T 

9pjj mVj 9P]j mVj 9pj^ mVj 9pj^ 

9M _ ^*^^m QSc ^ ^^Y,T 

9Pk ly 9Pk 9Pk ly 9Pk 

Note that 



V. 


Pk = q 
Pk^q 


(21a) 

(21b) 

(21c) 


Substituting equation (21a) into each of equations (18) in sequence results in 


Xg=- 


qjS 

mVj 


1 

,96 y/j mVi\ 96 J 


1 


(22a) 




qjS 

mVj 



Vl 

qi sinoi 


9£a 

9 n 



Vi /9 Fx,t\ 

-qi 


(22b) 


XA = 


qiSc 

2mVi^ 


?9Ca\ 

2Vi2 / 


. 1 

/9Fx,t\ 

cos a.l 

'8Fx^t\ 

U 

goC 1' 

/j 

mVj 

1, 

So 1\ 

/l 


(22c) 
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qjS g . ^ /8 Ca\ g sin 0t cos 0/8Fx,T 

^ ™ — i ; 5Yr-\ - ^ 


(22d) 


Xq = 


q^Sc 


2mVi" 


8<^ 

qc 


,9 


2V/ 


2Vi2 /9Ca 

— r— COS Q!t — 

g(jC 1\ 8n 


mVi 


'8F 


X,T\ 


/8F- 


8q 


+ — i cos O'! 
'1 go 


X,T' 


8n 


(22e) 


Xq = 


Xa = 


QlS 

BCj 

mVi' 

[all 

qiS 

Xl 

mVj 

So 


, 1 


qi 


COS 




(22f) 


mVi 


8F 


X,T 


ii + “A“5ry, ^ “A“»^ 


8F 


X,T 


+ Ri 


8F' 


X,T 


+ 2q. 


8F 


X,T 


8q 


(22g) 


q,S 

Xa = jj^sina 


dCA 1 . /3Fx^t\ 


/8Ca^ 

\ +-J— 1 


l»T 1 

'/mVi' 

^ 9T / 


(22h) 

(22i) 


8a \ 


_ q£_ ) 

h ~ m Vi ^ ai ^ah 8M 


Rl/8p\ Rl/8p' 

Pi 




8R 


1 99 


^8hj^^A, 


,>+_!_ A MiH f!!si, 

mVj \ v^/j\ 


Rl/8p\ Rl/8p 


9Fx,T 


<hjdp\ /9Fx,t 


Wh/^V 8q 


(22j) 


Equations (22) express the dimensional X derivatives in terms of the partial derivatives 
of with respect to the physical variables. Equivalent expressions in terms of the 
complete coefficient form of the stability derivatives are 
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c. + 


mVi -^5 mViV 96 


'9F 


X,T 



p g sin 6>i cos 4>(9 Fx^t\ 
^ 'mVj ^e~ goinVi \ an / 



+ Rj^ 



+ 2Qj 



1 




i = -M. 

mV, 


^9F- 


'AS 


mgo 


Sin a- 


X,T 


an 


Xt=-;2^Ca_+ * 


/9F- 


X,T 


mVi at mVi\ 9T 



(23a) 

(23b) 

(23c) 

(23d) 

(23e) 

(23f) 

(23g) 

(23h) 

(23i) 


(23j) 
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The expressions for the Z stability derivatives can be written simply by replacing 
X with Z, with Cjq’, and with F2^T “ equations (22) and (23). Simi- 

larly, expressions for the M stability derivatives are obtained by replacing X with 
M, with -cCm, ^x,T ^Y,T> mVj with ly in equations (22) and (23). 

The relationship of the stability derivatives developed herein to those discussed in 
references 2 and 3 is of interest. The derivatives with respect to physical variables as 
used in this paper are, in many cases, analogous to, and may be equal to, the derivatives 
termed "zero mass derivatives" in references 2 and 3. The stability derivatives of this 
paper (i.e., derivatives with respect to the motion variables) are analogous to those 
termed "equivalent elastic derivatives" in references 2 and 3. The values of the equiv- 
alent elastic derivatives as defined in references 2 and 3, however, will be different from 
those of the stability derivatives defined herein because the equivalent elastic derivatives 
were defined under the assumption of constant dynamic pressure. This assumption 
results not only in the omission of the dynamic -pressure contribution to the u deriva- 
tives, but also in a completely different form for the normal -acceleration contribution 
to all derivatives. Note in particular that for an vinaccelerated reference flight condition 
(qj = 0^, equations (18) show that the derivatives with respect to angle of attack reflect 
no contributions from normal acceleration. In contrast, all the equivalent elastic deriva- 
tives of references 2 and 3 are influenced by normal -acceleration contributions termed 
"inertial relief." 


Formulation of Equations of Motion 

The equations of motion are formed by expressing the conditions for dynamic equi- 
librium in each degree of freedom (eqs. (1)) in terms of the motion variables. The sta- 
bility derivatives defined in the preceding section represent the aerodynamic and thrust 
contributions in terms of the motion variables. The inertial contributions must now be 
expressed in terms of the same motion variables. For the airplane as a whole, the iner- 
tial reaction to gravity and translational accelerations in vector form is 



where 


g = i (-g sin e) + j (g cos e sin 0) + ?(g cos 0 cos (j>) 
V = iu + jv + k^w 
^ = Tu + Jv + kw + w X V 
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and 


cD = Tp + Tq + 


and i , j , and k are unit vectors directed along the X, Y, and Z body axes, 
respectively. 

Similarly, the inertial reaction in the airplane rotational degrees of freedom is 

Mj= -THx - iHy - kHz - w X H (25) 

where the angular momentiim H can be written H = I • w and the components of the 
inertia dyadic I in the body -axis system can be written in matrix form as 


1 = 


ix 

0 


0 

lY 


-I- 


xz 


-Ixz 

0 


if inertial symmetry about the X-Z plane is assumed. 

After performing the indicated vector operations, the inertial components in the 
six airplane degrees of freedom are found to be 


FXji = m(-g sin 0 - u - qw + rv) (26a) 

Fy^I = m(g cos 0 sin <p - v - vu + pw) (26b) 

^Z,I = ’^(S ® 0 - w - pv + qu) (26c) 

^X,I “ ~^XP ^XZ^ ■ (^Z ■ ^y)*!^ ^XZ^P (26d) 

My,i = -lytl - (ix - Iz)p^ - ^Xz(P^ - (26e) 

^Z,I == -^Z^ ^XZP " (■’•Y “ ^x)p‘l “ ^XZ*!^ (26f) 


In order to relate the variables in equations (26) to the reference and perturbation 
values of the motion variables adopted in the preceding section, use will be made of the 
defining equations for angle of attack and sideslip, 
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(27a) 


a - tan"l ^ 

/3 = sin-1 V (27b) 

the expression for altitude rate, 

h = u sin 0 - V cos 0 sin (^ - w cos 0 cos <p (28) 

and the expressions for the Exiler angle rates, 

0 = q cos (p - r sin p (29a) 

<J) = p + (q sin p + r cos p)taxi 0 (29b) 

p = {q, sin p + r cos p)sec 0 (29c) 


Note that equations (27) relate the linear velocity components to aerodynamic flow 
angles and equations (29) relate the angular velocity components to an earth -fixed refer- 
ence. If a quiescent. atmosphere and a nonrotating earth are assumed, the linear and 
angular velocities of equations (27) and (29) are identical to the velocities relative to 
inertial space appropriate for use in equations (24) and (25). 

The acceleration sensed by a seismic translational accelerometer located at the 
airplane center of gravity is Fi/m. The Z -component of this acceleration expressed in 
g units is defined herein as the normal acceleration n. Thus, from equation (26c) 

n = ^ cos 0 cos p + ^(-w - pv + qu) (30) 

&o ®o 

It should be noted that n is directed along the Z-axis and is, therefore, not iden- 
tical with the load factor, which is usually defined as the component of inertial force nor- 
mal to the flight path and the Y-axis, relative to the airplane weight, 

F<7 I cos a - Fv I sin a 

Load factor = — — — — (31) 

in go 

Equations (24) to (31) are applicable to large amplitude airplane motions in six 
degrees of freedom but are, in general, nonlinear. The linearized equations in the three 
longitudinal degrees of freedom usually have been derived by constraining the asymmetric 
motion variables (v,p,r,0) and their time derivatives to zero and expressing the remaining 
variables as the sum of a steady reference value and a small time -dependent perturbation. 
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The equations are then linearized by neglecting squares and products of perturbations. 

On examining the impact of the above constraints on equations (29a) and (26a and c), it is 
concluded that the only reference flight conditions which are both symmetric and steady 
are those with zero pitching velocity, implying a straight flight path. If an altitude - 
dependent atmosphere is considered, the admissible reference flight conditions are fur- 
ther constrained to straight level flight. 

It is apparent that analysis of flight at load factors significantly different from unity 
requires either violation of the small perturbation assumption or relaxation of either the 
symmetric constraint or the steady constraint on the reference flight condition. In order 
to retain small perturbations for elevated load factors, a symmetric, pseudosteady ref- 
erence flight condition having constant pitching velocity has been assumed on occasion. 

This condition is visualized as a wings -level pullup with a vertically curved flight path. 
Reference values of 0, h, and h can be defined only at one instant of time. Physical 
interpretation of the perturbations at other instants then becomes unclear under this 
assumption. 

Alternatively, a steady, pseudosymmetric reference flight condition can be assumed. 
Steady flight is achieved by constraining the time derivatives of all variables to zero, 
which requires that the angular velocity vector cS be directed along the gravity vector g. 
This condition can be visualized as a constant -altitude banked turn, or if a uniform atmo- 
sphere is assumed, a climbing or diving spiral. Approximate symmetry is achieved by 
constraining only the sideslip angle to zero. The reference flight condition must then 
satisfy the conditions for static equilibrium in all degrees of freedom, and the linearized 
longitudinal equations of motion are formulated by allowing small perturbations in only 
the longitudinal degrees of freedom. The equations of motion developed herein will be 
formulated for this steady, pseudosymmetric class of reference conditions since the truly 
symmetric conditions are special cases of this class. 

By applying equations (27) to the case of zero sideslip, the velocity components can 
be written as 

u = V cos a 

V = 0 

w = V sin a 

The longitudinal variables are stated in perturbation notation as 

V = Vi(l + u) (See eq. (17).) 
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a = ai + 0!p 

0 = 01 + 0p 

q = qi + qp 
h = hi + r* h dt 

h = hi + hp 

Then, for small perturbations 

u = Vi cos O’! + uVi cos «! - «pVi sin 0 !i 
w = Vi sin ai + uVi sin a-^+ ofpVi cos a-^ 
u = uVi cos «! - aVi sin ai 
w = uVi sin «i + aVi cos 
and, by substitution in equations (26), (28), and (30), 
Reference value 


Perturbation 


Fx,i = m^-g sin 0i - ViQj sin ai^ + mVi^-uqi sin 0 !i - u cos o!j - otpqi cos otj 

— ' (32a) 




+ a 


bi sin a 1 - Qp sin a i - 6p ^ cos 0i 


Fy^I = ni^g cos 01 sin tp - Vir cos + mVi sin o!i - r cos aij + «p^P cos 


+ Vip sin CKi^ 


+ r sin 0 ! 1 ^ - 0p ^ sin 0i sin <p 


(32b) 


FZ, 1 = m^g cos 01 cos (j) + Viqi cos a + mVi^i cos «! - u sin - otpqi sin a ^ 

- a cos «i +qp cos «i - 0p ^ sin 0i cos ^ 


(32c) 
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Reference value 

Perturbation 


r 

^X,I = 'pz ■ + ^XZP^l 

\ f 

+ <ip |;(iz ■ ^y)^ + ^xzp] 

(32d) 

My,I = -(^X - ^z)pr - Ixz(p^ - 

-qly 

(32e) 

Mz,I = -(ly - Ix)p^i - Ixz*!!^ 

+ qp [(ly - Ix)p - Ixzr] 

(32f) 


h = Vj^cos otj sin - sin cos cos 0) + uVi^cos sin - sin aj cos 0j cos 0^ 

- Q!pV]^(sin o!j sin + cos 0 !j cos 0j cos 0^ 
+ 0pVj^cos Q!j cos 0j + sin sin 0^ cos 0^ 


g Vi 

n = ^ cos 01 cos 0 + Ql cos Q! i 

On on '*• 


(33) 

ViL X . 

+ cos Q!j - u Sin Q!j - QfpQi sin cnj 

- a cos Q!]i + Qp cos aj) - dp ^ sin 0j cos 0 


(34) 


Unfortunately, equations (32b, d, and f) show that in turning flight, the lateral - 
directional components of the inertial reactions depend on the longitudinal motion pertur- 
bations. Clearly, the resulting lateral -directional perturbations could couple back into 
the longitudinal motions. In order to perform a longitudinal analysis without detailed 
knowledge of lateral -directional characteristics, the lateral -directional excitation will be 
neglected, with the knowledge that the error introduced diminishes to zero as the turn 
rate is reduced to zero. 

For the steady, pseudosymmetric reference condition, assume that static equilib- 
rimn of rolling and yawing moments can be achieved at zero sideslip with no net lateral 
force arising from aerodynamic and thrust sources. Equation (32b) then becomes 

^ cos 0j sin 0 - r cos Qfx + P sin = 0 (35) 

and from equations (29a and b) 


qj cos 0 - r sin 0 = 0 

(36a) 

p + ^qj sin 0 + r cos 0^ tan 0^ = 0 

(36b) 
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static equilibrium in the longitudinal degrees of freedom is expressed by substitution of 
the reference portion of equations (32a, c, and e) into equations (2) 

x(6i,Q!i,qi,ni,Mi,Ri,qi,Ti) = :^ sin 0^ + qj sin aj (37a) 

z(6i,Q!i,qi,ni,Mi,Ri,qj,Ti) cos cos 0 - q^ cos orj (37b) 

M(6i,ai,qi,ni,Mi,Ri,qi,Ti) = pr + - ^2) (37c) 

The system of eight equations made up of equations (35) to (37) and the reference portion 
of equations (33) and (34) is sufficient to define 6j, otj, q^, Tj^, 6-^, 0, p, and r 
if values of Vj, Mj, Rj, q^, n^, and hj are specified. It is interesting to note 
from equation (37c) that an inertial pitching moment exists in steady turning flight, which 
can influence the control-limited load factor of interest for fighter airplanes. 

The system of linearized longitudinal perturbation equations can now be formed 
from equations (2) by expressing the aerodynamic and thrust contributions in terms of the 
dimensional stability derivatives, expressing the inertial contribution by equations (32a, 
c, and e), supplementing the system with the altitude equation (33), and subtracting the 
reference equations (37) and the reference part of equation (33). The resvilting equations 
are 


- cos + u^X(i - qj sin + sin + «p(Xq, - qj cos oijj 

+ qX^ + qp(Xq - sin Oi) + 0p^X0 - cos + hpXh + SpXs + TpXx = 0 (38a) 

u^Z^ - sin + u^Z0 + qj cos ajj + - cos + a^i^a ~ CLj sin oijj 

+ + Qp^Zq + cos + 0p^Z0 - sin cos 0^ + hpZjj + 6pZg + TpZ>p = 0 (38b) 

uMg + uM^ + (xM^ + ttpM^ + q^M^ ~ *lp^ ^P^S ^p^h + 6p^6 + TpMT = 0 (38c) 


uVj^^cos Oj sin - sin <Xj cos cos 0j - tXpVj^sin sin + cos cxj cos 0^ cos 

+ 0pVifcos oj cos 01 + sin oii sin 0j cos 0^ - fip = 0 


(38d) 
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To obtain a solution, this set must be supplemented by equations for the two control sys- 
tems, which will not be stated here, and by the perturbation equation obtained by subtract- 
ing equation (36a) from equation (29a), that is, 

0p - qp cos (j) = 0 (38e) 


since = 0. 

For the special case of a steady, straight, symmetric reference flight condition, 

(^ = p = qj^ = r = 0 and the equations of motion are simplified accordingly. In particular, 
equation (38d) becomes 

uVj sin^0j - - CKpVj cos^0j^ - + 0pVj cos^0j - ofi) - hp = 0 (39) 


and the reference equations are 


Xl = ^ sin 01 
Zl= -^cos 01 


Ml = 0 \ 

hi = Vi sin^0i - ai^ 

ni = ^ cos 0-1 
&o 


(40) 


If the reference condition is further constrained to level flight, then 0i = «i, and the 
resulting simplification of the equations of motion is obvious. 

It should be noted that equations (38) can be considered to be a hybrid set of equa- 
tions of motion in the sense that the motion variables are resolved on stability axes (uVi 
is directed along Vi), whereas the X and Z forces are resolved on a body -fixed axis 
system whose orientation in the airplane plane of symmetry is arbitrary for rigid air- 
planes and is determined by structural consideration for elastic airplanes. The more 
conventional stability -axis system of equations of motion is obtained from equations (38) 
by simply setting ori = 0. K the stability -axis system is to be used for elastic airplanes, 
it is recommended that the stability derivatives be determined first on a structurally 
appropriate body-axis system and then be transformed to the stability -axis system by 
the usual force transformation. 
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Calculation of Stability Derivatives 

In this section, an influence -coefficient procedure for calculating the stability 
derivatives is outlined. This procedure has been implemented in a digital computer pro- 
gram which is presented and described in the appendix. In addition to the assumptions 
already made in defining the stability derivatives, the following assumptions are made in 
order to simplify the calculation procedure: 

1. Aerodynamic loadings are derived from a linearized, inviscid, thin lifting -surface 
theory. 

2. Singular aerodynamic loads arising from edge suction are neglected. 

3. All unsteady aerodynamic effects are neglected. 

4. All contributions of the propulsion system are neglected. 

The effect of assumptions 1, 3, and 4 is to eliminate the variables R, a, and T 
from the list of physical variables in table I. Equations (6b) and (7b) are then represented 
by the matrix expression 

{^) = [A] (c; + e> (41) 

where is an aerodynamic influence -coefficient matrix which can be obtained from 

any suitable lifting -surface theory for a given Mach number and paneling geometry. ' In 
the present development, aerodynamic symmetry is assumed, and each element of (a) 
represents the aerodynamic normal force acting at a load point (x£,yf,0) due to a unit 
change in mean -camber -surface slope at a symmetric pair of control points ^Xa,±ya»0) 
for unit dynamic pressure. 

The structviral operator ^ in equation (10) can also be represented by an influence- 
coefficient matrix such that 



where each element of [S] is the elastic change in mean-camber-surface slope at an 
aerodynamic control point (xa,yaj0) due to a unit force applied at a symmetric pair of load 
points (xf,±yf,0). Combining equations (8a), (9), (11), (41), and (42) and solving for the 
aerodynamic -force set give 


2Sj 
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The coefficients Cn and Cm can now be written 



(44a) 

(44b) 


Under the second assumption stated at the beginning of this section, aerodynamic 
edge suction is neglected, and the aerodynamic force on each panel then acts normal to 
the mean camber surface at that panel. If the mean -camber -surface slope is small 
everywhere, the axial -force coefficient is approximately 



(44c) 


where ^ is the set of mean-camber -surface slopes at the panel load points 

® ' •&«,() 

It should be observed, however, that the neglect of edge suction implies a rather gross 
oversimplification of the true aerodynamic loads, particularly in subsonic flight, and the 
resulting axial -force coefficient must be considered incomplete. Despite this shortcom- 
ing, the axial stability derivatives will be formulated herein to provide a basis for future 
refinement and to allow an examination of their possible significance to the stability of 
elastic airplanes. 

The stability derivatives are calculated by substituting the partial derivatives of 
equations (44) with respect to each of the physical variables a, 6, qc|2Vi, n, q, M, 
and q into the coefficient form of equations (18). By denoting any physical variable 
by P, 
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(46a) 



The right-hand side of equation (43) can be interpreted as the product of an aero- 
elastic correction matrix, the aerodynamic influence -coefficient matrix, and a vector of 
pseudorigid aerodynamic boundary conditions which includes structural distortion from 
only the inertial loads. To simplify the notation, the aeroelastic correction matrix can 
be defined as 


[l] -qC^[s] 


n-1 


(47) 


and the vector of boundary conditions as 

= {'^ig> + « (l) + 6 {^6} - |ngo {1} + (q - Pr) [m] (xf^ 

(48) 

At the reference flight condition, 

(‘^>1 = ("jig) + «l(l> + + &:^lgo [m] {1} - pr [ni| {xj^ 

(49) 

= (50) 

and 


{"f> 1 = ^ (^6,0 Kl ® 1 1 1 ""1^0 ^ (1) - pr 


(51) 
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The partial derivatives appearing on the right-hand side of equations (46) can be 
evaluated in a straightforward manner with respect to the first five physical variables a, 
6, qc/2Vj, n, and 4 . A finite -difference technique has been chosen to evaluate the 
derivatives with respect to M and q. For increments in M and q which are small 
relative to the reference -condition values, the following matrices are defined: 

[[ 1 ] - (qj + 4q)[A]i[S]j‘ - [H - (qj - Aq)[AJi[s]j' 

N “ m 




H - qi|^M,l+AM]rs]J |Am,1+Am] - [H ■ ^l[AM,l-AM][S]j [^m,1-Am] 


2 AM 


The matrices of the partial derivatives of F2 ^a/^ ^nd ej with respect to each phys- 
ical variable are expressed in table in. These values along with the reference -condition 
values of equations (49) to (51) are used in equations (46) to find the derivatives of Caj 
Cn, and Cjji which are in turn substituted into the coefficient form of equations (18) to 
find the stability derivatives. 

fdFz aJ^'] 

TABLE m.- THE PARTIAL-DERIVATIVE MATRICES ( §p^ / 


Physical 
variable, p 

|sI'z,A/9j 

ipp. 


a 

01 [A] 1 ( 1 ) . 

«i[S(]Cb]iCA]i(1> 

5 



qc 

2Vi 

[B]i[A]i(^) 

«lC%]®llA]i{-^) 

n 

got®] i[s] [m] { 1 } 

«oKlM ^4irB]iCA],[s][m] (1) 

• 

q 

(xf) 

i 

[Sf] p.] + q 1 [b] J [A] J [s] [m] (xf) 

M 

i 

|^mAm]^)i 

^i(^|?mAm]{c^}i 

q 

1 

[B^[A],{cr>i 

S1[9i[B9] +01 

C^iWi 
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Evaluation of equations (49) to (51) requires knowledge of the mean-camber -surface 
slopes of the jig shape ^ as well as the reference -condition values 

of Sj, qj, nj, p, and r. Assume that the airplane configuration has been fully 

defined at a single, straight, level design flight condition and that the design air -load dis-< 
tribution is known. By applying equations (9), (42), and (45), the jig shape slopes are 
foimd to be 


{j^ig} “ {^es} “ ®des{56} " C^^des|^ ^ ^esSo ^de^ '{l}. 


^jig,f} “ {^es,f} " ®des{%f} 


+ "desSo 


“des](l> 


The values of the physical variables at any desired reference flight condition can be foimd 
by applying equations (33) to (37). Some simplification results, however, from the present 
assumptions under which the effects of thrust and nonlinear aerodynamics are neglected. 

If the thrust is assumed just sufficient to balance the X forces, equation (37a) can be 
neglected. Equations (37b and c) can be restated as 


^N,jig+“1 Qa as 2Vi /qc\"*' 1 8n “ 3q 


mVi / g ^ 

-=-TT ^ cos 01 cos (i) + qi cos 0 !i 

qiS \yi 1 1 1 


'in, Jig 




^^m . *11*^ ^^m . .. ^^m .... ^^m 


= sfe (iz - - Ixz(p2 - r- 


where 


jig “ “ |{0 [Ml jK 1 {ejig} 
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and 


Cm, jig = - ^ [b] 1 K 1 {€jig) 


The derivatives with respect to q are used to express the pr contributions since the 
inertial -force distribution for unit q is shown by equation (11) to be equal to that for a 
xmit value of -pr. For specified values of Mj and the derivatives appearing in 
equations (53) can be calculated by using equations (46) and table IH. The reference flight 
condition can then be fully defined by solution of equations (33) to (36) and (53). 

Effect of Structural Fixity Point Location 

The influence -coefficient procedure for calculating stability derivatives of an elas- 
tic airplane, as outlined in the preceding section, makes use of a structural influence - 
coefficient matrix each element of which is the elastic change in mean-camber - 
surface slope at an aerodynamic control point due to a pair of unit forces applied at a 
symmetric pair of load points. Evaluation of this matrix requires that the structure be 
assumed clamped at a structural fixity point. Any nonequilibrium distribution of loads 
applied to the structure must then be reacted by a force and moment imposed at the fixity 
point, and the deformed shape, calculated by using the matrix [s], will include implicitly 
the effects of these reactions. The deformed shape depends, therefore, on the fixity point 
location. Since the load distribution arising from a unit pertimbation of any single phys- 
ical variable is, in general, not in equilibrium, the values of the stability derivatives cal- 
culated by the procedure outlined will depend on fixity point location. 

When a complete and consistent set of stability derivatives is used in the equations 
of motion derived from the conditions for dynamic equilibrium, the reactions at the fixity 
point vanish and the deformed shape at any instant is, in a coordinate -free sense, inde- 
pendent of fixity point location. Some of the motion variables, however, for example, a, 

9, and their time derivatives, are not coordinate -free but are defined relative to the ref- 
erence axis system. Moreover the elastic surface slope changes defined by the matrix 
[S] are also expressed relative to the reference axis system and must be zero at the 
fixity point. Thus, the orientation of the reference axis system must be considered 
invariant with respect to a material element located at the fixity point, as indicated in 
figxire 1. 

It is apparent, then, that after application of the equations of motion, any quantities 
that are independent of reference axis system orientation in the airplane, such as char- 
acteristic roots, the variables u, 6, h, and the magnitudes of resultant force, moment, 
or acceleration vectors, are also independent of the fixity point location. Those quan- 
tities dependent on reference axis system orientation, such as the variables a, 9, q. 
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and components of force or acceleration vectors, are also dependent on the location of the 
fixity point assumed in the definition of the matrix 

ILLUSTRATIVE ANALYSIS 
The Example Airplane 

To illustrate the calculation and use of the stability derivatives defined in this 
paper, a longitudinal quasi-steady aeroelastic analysis of a high-performance supersonic 
transport airplane configuration has been performed. Although the analysis was per- 
formed in U.S. Customary Units, the results are presented in the International System of 
units (SI). The configuration selected was that used in reference 5 to illustrate a number 
of aerodynamic design integration features. This configuration was also the subject of 
a preliminary design study by The Boeing Company during which airplane weight and bal- 
ance information and elastic properties of the structure were developed. The author 
gratefully acknowledges the cooperation of The Boeing Company in supplying a symmetric 
structural slope influence -coefficient matrix and weight distribution data for use in the 
present analysis. It should be emphasized that these data resulted from a preliminary 
design study and are not necessarily applicable to an airplane design meeting all the 
requirements of a commercial supersonic transport airplane. In particular, the struc- 
tural elastic characteristics do not reflect any additional stiffening which might be neces- 
sary to achieve adequate controllability and flutter margin. 

For the present study, the aerodynamic effects of the vertical tails and engine 
nacelles were neglected and a new mean camber svu*face was designed to minimize the 
drag due to lift and to trim the airplane with zero control -surface deflection at M = 2.7, 
Cn = 0.07914, and a center -of -gravity location at 0.417c. The configuration planform 
and the paneling geometry used for subsonic and supersonic analyses are illustrated in 
figure 2. Corresponding aerodynamic influence -coefficient matrices were calculated by 
using the vortex -lattice method of reference 6 for the subsonic analysis (m^ = O.s) and a 
procedure based on the method of reference 7 for the supersonic analysis (Mj * 2.7). 

The load points and aerodynamic control points were located at the quarter -chord and 
three -quarter -chord points, respectively, of the panel midspan chords for the subsonic 
case and at the panel centroids and 0.95 panel chord points, respectively, for the super- 
sonic case. 

The structural influence -coefficient matrix supplied by The Boeing Company was 
transformed to the appropriate forms for the subsonic and supersonic panel arrangements 
by using a second-order, bidirectional interpolation procedure applied successively to the 
slope -point and load -point transformations. The fixity point of the structural influence - 
coefficient matrix used in the present analysis was located at 0.33c. 
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The weight distribution supplied by The Boeing Company was transformed by dis- 
tributing each original weight element among the three closest load points of the new 
panel layout in such a way that the moments about the X- and Y-axes were imchanged. 

The distribution of structure, contents, and payload weight, totaling 1.56 MN (351 560 lb), 
was considered invariant. Fuel capacity constraints, based on the fuel weight distribu- 
tion supplied by The Boeing Company, were assigned to each panel. Within these con- 
straints, fuel weight distributions yielding the following combinations of airplane weight, 
center -of -gravity location, and pitching moment of inertia were developed: 


Weight 

( 1 . 

ly X 10“6 

MN 

ib 

kg-m2 

slug -ft 2 

2.67 

600 000 

0.41 

61.61 

45.44 

2.67 

600 000 

.47 

62.43 

46.05 

1.87 

420 000 

.41 

56.33 

41.55 

1.87 

420 000 

.47 

57.43 

42.36 


The weight distribution for the design condition, W = 2. 10 MN (472 500 lb) and center of 
gravity at 0.417c, was also developed. 

The chordwise distributions of mean -camber -surface slope for the design shape 
e^jgg and the jig shape ejjg at several spanwise stations are shown in figure 3. These 
slopes are referred to a reference axis system oriented in such a manner that the angle 
of attack at the design condition is 2.8°. This orientation was chosen so that the average 
slope of the design mean camber surface was approximately zero. 

The longitudinal characteristics of the example airplane have been analyzed for a 
series of reference flight conditions including center -of -gravity locations of 0.41c and 
0.47c, the two airplane weights noted previously, and either five values of dynamic pres- 
sure ranging from 7.2 kN/m^ (150 Ib/ft^) to 47.9 kN/m^ (1000 Ib/ft^) for M = 0.8 or 
four values of dynamic pressure ranging from 14.4 kN/m^ (300 Ib/ft^) to 47.9 kN/m^ 
(1000 Ib/ft^) for M = 2.7. All reference flight conditions corresponded to straight level 
flight. For comparison, the rigid airplane having the design mean camber surface was 
also analyzed. 

Longitudinal control and trim were achieved by equal deflection of the inboard and 
outboard trailing-edge control surfaces outlined in figure 2. Because the control surface 
edges were not coincident with panel boimdaries, control deflections were represented by 
assigning to each panel a value of eg equal to the ratio of the control surface area lying 
within that panel to the total panel area. For each panel, the values of eg f and eg 
were assumed to be identical. 
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I III 


Aerodynamic Derivatives 

The computer program described in the appendix was used to calculate the longi- 
tudinal stability derivatives together with parameters describing the reference flight 
conditions. The values of angle of attack, control deflection, and normal- and axial -force 
coefficients calculated for the reference flight conditions are plotted in figure 4 for the 
rigid and elastic airplanes at the two center -of -gravity locations. The axial -force coef- 
ficient includes an increment of 0.00550 to represent skin friction and form drag at 
M = 0,8 or an increment of 0.00533 to represent skin friction and wave drag due to vol- 
ume at M = 2.7. The aeroelastic effects on trim angle of attack, control deflection, and 
axial -force coefficient are apparent from figure 4. 

The partial derivatives of Ca, Cn, and Cm with respect to the physical vari- 
ables are shown in figures 5, 6, and 7. Because unsteady aerodynamic effects were 
neglected, no partial derivatives with respect to a were calculated. When the stability 
derivatives are formed by using equations (18), however, stability derivatives with respect 
to dt will exist for the elastic airplane because of the contribution of the partial deriva- 
tives with respect to n. For the rigid airplane, all derivatives with respect to n, q, 
and q are zero. All Ca derivatives and the derivatives of Cjq- and with 

respect to Mach number depend on the reference values of a and 6 and consequently 
vary with dynamic pressure. Partial derivatives of Cjj and Cm with respect to the 
remaining physical variables are independent of the reference flight condition for the 
rigid airplane. 

The effects of aeroelasticity on the partial derivatives are apparent from figures 5, 
6, and 7. Note in particular the large destabilizing aeroelastic effect on dC-^jda. 

The Characteristic Stability Roots 

The longitudinal stability derivatives were calculated by applying equations (18) to 
the partial derivatives of figures 5, 6, and 7. The three-degree-of -freedom linearized 
equations of motion (eqs. (38)) were then solved for the characteristic roots and the trans- 
fer functions u/6p and np^6p. The perturbation part of equation (34) was included in 
the system of equations to define np. The altitude -dependent parameters required for 
the solution are shown in figure 8 and are based on the 1962 U.S. Standard Atmosphere 
(ref. 8). Solutions corresponding to a locally uniform atmosphere were also obtained by 
setting the derivatives X^, Z^, and Mjj equal to zero. 

The real and imaginary parts of the characteristic roots are presented as functions 
of dynamic pressure for the rigid airplane in figure 9 and for the elastic airplane in fig- 
ure 10. Real roots are shown as faired lines without symbols. The calculated real and 
imaginary parts of complex roots are indicated by symbols, with the same symbol shape 
used for the real and imaginary parts of any particular conjugate pair of complex roots. 
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Calculations were made at only four or five values of dynamic pressxire, and therefore, 
the lines faired between calculated points should be considered as only illustrative. In 
particular, the locations of branch points, that is, points of transition between a pair of 
real roots and a pair of complex roots, are only approximate. 

The system of equations yields five characteristic roots including one root arising 
from the differentiation of h in the altitude equation (38d). For the cases representing 
a locally uniform atmosphere, the altitude equation is uncoupled from the remaining equa- 
tions and one root is always zero. With a standard atmosphere, the altitude root was 
usually too small to be apparent in figures 9 and 10 except for the aft c.g. cases at 
M = 0.8 (figs. 9(b) and 10(b)) where the altitude root can be seen at low dynamic pres- 
sure. The remaining roots in most cases can be clearly identified as a short -period 
pair and a phugoid pair although these root pairs were not always oscillatory. When the 
remaining roots are all complex, they can be recognized clearly as a phugoid pair and a 
short -period pair, with the phugoid pair having the smaller magnitude (smaller undamped 
natural frequency). Identification of real roots is not as straightforward, but in most 
cases, the faired lines in figures 9 and 10 can be used to relate real roots to either the 
phugoid or the short -period mode at other values of dynamic pressure. 

The assumption of a locally uniform atmosphere had a destabilizing effect on the 
phugoid roots, usually manifested as a reduction in phugoid frequency, but had little effect 
on the short-period mode, particularly when this mode was oscillatory. A comparison of 
figures 9 and 10 shows that the effect of aeroelasticity on the short -period root pair was 
to reduce both frequency and damping with the frequency reduction being more pronoimced. 
As a result, the short -period mode of the elastic airplane in many cases became over- 
critically damped and therefore yielded two aperiodic (real) roots. This tendency was 
most apparent with heavy weight, aft c.g., high dynamic pressure, and subsonic speed. 

In all cases, however, the short-period root pair remained stable as indicated by the neg- 
ative real parts. 

The phugoid roots of the rigid airplane (fig. 9) were always small in magnitude and 
corresponded to either a lightly damped low-frequency oscillation or a mild aperiodic 
divergence. In contrast, the elastic airplane (fig. 10) exhibited, in many cases, phugoid 
roots which were much larger in magnitude and were strongly divergent in either an 
oscillatory or aperiodic manner. Again this tendency was most apparent with heavy 
weight, aft c.g., high dynamic pressure, and subsonic speed and appeared to be particu- 
larly sensitive to changes in weight or possibly weight distribution. 

In an attempt to identify the reasons for the aperiodic short -period roots and the 
strong phugoid divergence, the equations of motion were simplified to two second-order 
sets yielding decoupled approximations to the short -period and phugoid modes. For a 
locally imiform atmosphere and a wings-level reference flight condition, equation (38d) 


can be omitted and equation (38e) becomes dp = qp. After resolving the force equations 
onto stability axes and applying the Laplace transformation, equations (38a, b, and c) can 
be written in matrix form as 


(1 - xa^s) - Xa,s 

■®^Q!,S " ^Q!, 

“®^u,s “ 2^u,s 

s(l - Z^^s) - Z, 

-sMa - M(i 

-sM^ - 


-s^Zq^s " + ^q,s) “ ^0,s 

s2(l - M-) - sMq - M0 


u 


«p) 


Xs,s 

{Z6,s 


, 


(54) 


A decoupled approximation to the short -period mode is obtained by constraining u to 
zero and assuming Zg g and Mg to be negligible 


s(l - Z^^s) - Zqj^s -sZq^g - ^1 + Zq^s) 

-sM^ - Mq, s(l - M^) - Mq 

The decoupled short -period roots are then the roots of 

s2 + 2^gpO;gpS + ^sp^ “ ® 

where 

2_ Z«,sMq - Mc,(l + Zq,s) 

“(l-Z4s)(l-Mq)-Z^^sMi 

_ (1 - Zajs)Mq - (l - Mq)Zo;,s " Zq^sMa - Mi(l + Zq^s) 

(1 - Z4,,)(l - M^) - 



(55) 


(56a) 


(56b) 


A decoupled approximation to the phugoid mode is obtained by retaining only the zero- 
order terms in the moment equation and assuming Xq^s and Zq^g to be negligible 
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The decoupled phugoid roots are the roots of 

s2 + 2?pWpS + Wp2 = 0 

where ojp^ and 2?pa>p are found by expanding the determinant of the matrix in 
equation (57). 

Figure 11 presents a comparison of the decoupled root parameters from equa- 
tions (55) and (57) with those from the coupled system (eq. (54)) for the elastic airplane 
with forward center -of -gravity location at M = 0.8. The coupled and decoupled short- 
period parameters are in good agreement except for those cases in which the short-period 
mode was aperiodic. The decoupled phugoid approximation, however, was generally poor 
and, in particular, failed to yield the large phugoid magnitude shown by the coupled solu- 
tion for the higher dynamic pressures with heavy weight. Observe, however, that the 
total damping 2^spWsp + 2^pCOp was well represented by the decoupled solutions at high 
dynamic pressure. This fact implies that the coupled system provides a mechanism for 
transfer of damping from the phugoid into the short -period mode. Examination of other 
simplifications of equation (54) has led to the conclusion that the root-coupling mechanism 
is provided primarily through the derivatives Mq and and becomes important when 
the numerator of equation (56a) - Mqj^I + Zq^s) is small. 

It should be noted that for the elastic airplane, the inertial contributions to the 
derivatives Mq, Zq, M^, and arising from normal -acceleration effects (eqs. (18)) 

were much larger than the aerodynamic contributions and were responsible for large neg- 
ative values of Mq and Zq and large positive values of Mq, and Zq,. As a result, 
the value of t 0 sp 2 remained positive (stable) for the lightweight cases and near zero for 
the heavyweight cases of figure 11 in spite of the large unstable values of M^^. The 
effect of weight on Wgp2 resulted primarily from the larger magnitude of Z^y arising 
from the reduced mass of the lightweight condition. 

Static -Stability Considerations 

Static longitudinal stability can be interpreted as the tendency for either speed or 
normal acceleration to be restored to its equilibrium value following a disturbance. The 
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two t 3 ^es of static stability to be discussed, therefore, will be termed speed stability and 
maneuvering stability. The concepts of static margin and maneuver margin have been 
used as direct quantitative measm-es of the two types of static stability. These stability 
margins are visualized as the distance from the c.g. to either the neutral point or the 
maneuver point where these points are the c.g. locations for which either the speed sta- 
bility or the maneuver stability is zero. The stability margins can be quantified, there- 
fore, as values of the total derivative dCm/dCiq- evaluated under constraints appropriate 
to a disturbance in either speed or normal acceleration. Approximate expressions for 
the stability margins in terms of the stability derivatives exist in the literature, but their 
applicability to elastic airplanes should be examined. 

Other measures of static stability express the observable consequences of the inter- 
action of static stability with other airplane properties. For example, the undamped nat- 
ural frequencies of the pht^oid and short -period modes reflect the interaction of the speed 
stability and maneuvering stability, respectively, with the inertial properties of the air- 
plane. Similarly, the static control parameters 6p/u and Sp^np are measures of 
speed stability and maneuvering stability, respectively, relative to control effectiveness. 

In the present analysis, the static control parameters are determined from transfer func- 
tions obtained from solution of the equations of motion and, along with the stability mar- 
gins, will be utilized in describing the static stability. 

The asymptote form of the frequency response amplitude of the u/6p and npy^6p 
transfer functions is illustrated in figure 12 for the elastic airplane in the heavy condition 
with c.g. at 0.41c and M = 0.8 for five values of q. This form of plot, commonly 
termed Bode plot or corner plot, is made up of straight-line segments connecting points 
at the pole and zero frequencies of the respective transfer functions. For the cases 
shown, the values of logj^Q of the phugoid frequency range from about -1.25 at the low- 
est dynamic pressure to about -0.5 at the highest dynamic pressure. Note that in the 
u/6p responses, a plateau exists at frequencies just below the phugoid frequency in which 
the amplitude u/6p is independent of frequency. For the uniform -atmosphere cases, 
this amplitude level is equal to the zero -frequency asymptote and can be used to obtain a 
static value of the static speed control parameter 6p/u. The value of 6p/u thus 
obtained corresponds to a constant -altitude constraint because even though the final alti- 
tude following a step control input is different, in general, from the initial altitude, the 
altitude change has no effect imder the assumption of a uniform atmosphere. The sign 
of 6p^u is determined from the phase angle of the asymptote form of the frequency 
response (which is either 0° or 180°) in this same frequency range. 

For the standard -atmosphere cases, the zero-frequency asymptote of u/Sp has 
little significance because of the loi^ time required for the aircraft motions to approach 
the final conditions. The value of 6p/u corresponding to frequencies just below the 
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phugoid frequency, on the other hand, is believed to represent a meaningful measure of 
the tendency within a time range measured in minutes for the speed to return to its trim 
value under conditions of no control action, that is, no altitude constraint. Thus, a dis- 
tinction is drawn between two speed stability indices, both expressed as values of 6p/u. 
One, determined from the standard -atmosphere solutions, is a measure of the imcon- 
strained speed restoring tendency, and the other, determined from the uniform -atmosphere 
solutions, is a measure of the control required to change speed at constant altitude. 

The static maneuvering control parameter 6p|np is traditionally determined from 
solution of a two -degree -of -freedom (decoupled short period) set of equations of motion. 
The np/6p frequency response curves identified as "constant speed" in figure 12 were 
obtained from solutions of equation (55) combined with an appropriate form of the normal- 
acceleration equation. The asymptote form of the frequency response of Opy^Sp obtained 
from the standard-atmosphere and uniform -atmosphere solutions of the complete three- 
degree-of -freedom equations of motion are also shown for comparison. For the lower 
three values of dynamic pressure, the levels of np^6p in the frequency range between 
the phugoid and short -period frequencies from the standard- and imiform -atmosphere 
solutions are identical with the zero -frequency asymptote of the constant-speed solutions. 
Note that at these three values of q, the short -period mode is oscillatory and its fre- 
quency is well separated from that of the phugoid mode. For the two higher values of q, 
the short -period root pair has branched into two real roots, and the frequency range of 
the np^6p plateau between the phugoid roots and the smaller root of the short -period 
pair is greatly diminished or nonexistent. The plateau level of np^6p from the three - 
degree-of -freedom solution is no longer identical to the zero -frequency asymptote of the 
two-degree-of -freedom solution. While the static maneuvering control parameter Spy^Hp 
obtained from the three -degree -of -freedom solution is probably more representative of 
the actual case than that from the constant -speed solution, this distinction may be of only 
academic interest when the large dynamic instability of the phugoid mode for these high- 
dynamic -pressure cases (fig. 10) is recalled. 

The values of 6py^u and Spy^np obtained from the appropriate three -degree -of - 
freedom transfer fimctions by the procedure just described are presented in figure 13 for 
the rigid airplane and in figure 14 for the elastic airplane. In a few cases, the pole -zero 
configurations of the transfer fimctions were such that the appropriate plateaus in the fre- 
quency response curves did not exist, and these points are omitted from figures 13 and 14. 
The assumption of a locally uniform atmosphere had little effect on the maneuver control 
parameter 6p^np but showed a destabilizing effect on the speed control parameter 6p^u, 
particularly at supersonic speed. This result implies that the unconstrained speed- 
restoring tendency, indicated by the standard-atmosphere results, is significantly more 
stable than would be inferred from measurements of the control gradient with speed at 
constant altitude, as indicated by the uniform -atmosphere results. 
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The static speed control parameter and maneuvering control parameter for the 
rigid airplane reflect the expected reduction in static stability with rearward center -of- 
gravity movement and increase in control sensitivity with increasing q. Comparison 
of figures 13 and 14 reveals several manifestations of aeroelasticity. At low dynamic 
pressure, an aeroelastic reduction in static stability is indicated by reduced magnitudes 
of both Sp^u and Sp^np. The subsequent increase in the magnitude of Sp^Up at high 
dynamic pressures is related to the aeroelastic loss in the control effectiveness deriva- 
tive Cmg. The aeroelastic loss in Cmg is not sufficient, however, to explain the large 
increase in Sp^^u at high dynamic pressures shown in figure 14, particularly when the 
large destabilizing aeroelastic effect on the partial derivative shown in fig- 

ure 7 is recalled. 

Further insight into the nature of aeroelastic effects on static stability may be 
gained by an examination of approximate forms relating static margin and maneuver 
margin to the stability derivatives. If the atmosphere is assumed to be locally \miform, 
the reference flight condition to be straight and level with aj small, and the perturbed 
flight condition to be steady trimmed flight along a curved path, then 

= %,1 + + ‘Ip ^ ^Nq + ^p^Ng 

Cm = 0 = QfpCm^ + uCm^ + qp ^mq + SpCj^g 


Also, 


Vl 

Thus, by substitution 

+ ^(2Cn, 1 + Cn^) + np(|^ Cn^ - ^ + 6pCNg = 0 (58a) 

^m = ‘^pCm^ + ^^m^ + *^p 2 ^ ^*^q ^ ® (58b) 

The stability margins can be expressed as values of the total derivative dCj^/dCj^ 
for the vmtrimmed perturbed state (Cm o) obtained by constraining 6p to zero. The 
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form of an additional constraint is used to distinguish between static margin and maneu- 
ver margin. For level flight at varying speed, the constraint np = 0 is imposed and 


Static margin = ' — ^ 




'm. 


u 


Wnp=6p=0 2Cn,i 


(59a) 


Similarly, for maneuvering flight at constant speed, u is constrained to zero and 
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(59b) 


Alternatively, by constraining either np or u to zero, equations (58) can be 
solved directly for the static control parameters 


T=2Cn,i 
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The approximations provided by equations (60) are plotted in figures 13 and 14 and, in 
nearly all cases, agree very well with the results from the transfer function solutions for 
the uniform -atmosphere cases. Thus, the mechanism responsible for the extreme values 
of 6p/u observed for the elastic airplane at the highest dynamic pressure must be 
embodied in equations (59) and (60). Equations (60) imply that the aeroelastic reduction 


in Ct 


should affect both 


in a similar fashion. 


'mg oiiuiuiA cuic^k uutii 5-pl VL and ®p/^p 
in passing that the combined control effectiveness parameter 
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found to be independent of center -of -gravity position even for the elastic airplane 


It might be noted 

Cm,- 


^Q! 


was 


The difference between the trends of Sp^u and Spy^np with dynamic pressure 
arises from the difference between the static margin and maneuver margin. These sta- 
bility margins evaluated from equations (59) are shown for the rigid airplane in figure 15 
and for the elastic airplane in figure 16. Values of which is indicative of the 

distance from the center of gravity to the aerodynamic center, are also shown to demon- 
strate the contribution of and Cm- to the static margin and the contribution of 
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Cjjq and Cjnq to the maneuver margin. Considering first the rigid airplane (fig. 15), 
the maneuver margin is shown to be more negative (stable) than particularly 

at the subsonic speed. The generally stabilizing contribution of Cmq to the maneuver 
margin is thoroughly documented in the literature. The contributions of the u deriva- 
tives to the static margin are not as well documented but are shown by figure 15 to be 
significant for the cases considered. For the rigid airplane, the u derivatives arise 
only from the effects of Mach number perturbation. Their contribution to the static mar- 
gin is generally stabilizing at M = 2.7 but destabilizing at M = 0.8. For the subsonic 
cases with the rearward c.g. location, the static margin is shown to be positive (unstable) 
in spite of the value of Cma/CNo, of -0.03. This result is compatible with the existence 
of a divergent root of the phugoid pair shown for the imiform-atmosphere cases in fig- 
ure 9(b). Figure 15 also shows that the charge in static margin with a change in c.g. loca- 
tion is significantly different from the change in (^) . This result arises primarily 

from the effect of The contribution of Cfjq to the maneuver margin was found to 

be much less significant. 

For the elastic airplane (fig. 16), the values of became highly positive 

at high dynamic pressure. The u derivatives now include the aeroelastic effects arising 
from dynamic -pressure perturbations (eqs. (18)) which were found to predominate over the 
Mach number effects. Similarly, the q derivatives include the aeroelastic effects from 
Up which were much larger than the direct aerodynamic effects of pitching velocity. The 
resulting contributions of Cm- to static margin and of Cmq to maneuver margin were 
always stabilizing. Moreover, the effects of and Cjjq were significant and always 

acted to reduce the contribution of Cm^y'Cjj^ to the stability miargins. In particular, 
for the highest dynamic pressure at M = 0.8, the values of CN^y^2CN^X were of the order 
of -1; the contribution of equation (59a) was therefore of little significance 

and the static margin became extremely stable, leading to the large values of 6p^u pre- 
viously noted. 

Figure 16 also shows that as dynamic pressure was increased, the stability margins 
became less sensitive to c.g. location, and the trend of static margin with c.g. location 
actually reversed at high dynamic pressures for M = 0.8. Furthermore, calculations 
made at intermediate c.g. locations showed that the variations of both static margin and 
maneuver margin with c.g. location became nonlinear in such a manner as to cast doubt 
on the existence of a neutral point Or maneuver point. 

Added confidence in the validity of equations (59) is gained by correlating the resiolts 
shown in figure 16 with the characteristic roots for the uniform -atmosphere cases of fig- 
tire 10. Note that a divergent real root is shown in figure 10 for each case for which 
either the static margin or maneuver margin is unstable as shown in figure 16. Further- 
more, the appearance of an aperiodic short -period mode at intermediate values of dynamic 
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pressure for the lightweight, rearward c.g., subsonic condition (fig. 10(b)) appears to be 
related to the minimum in the stable values of maneuver margin shown in figure 16(a) at 
intermediate values of dynamic pressure for the corresponding condition. This latter 
correlation is suggested by the similarity between the form of equation (59b) for maneu- 
ver margin and that of equation (56a) for Wsp^. 

Effect of Dynamic-Pressure Derivatives 

One of the departures from past practice embodied in the present development is 
the inclusion of the aeroelastic effects of dynamic -pressure perturbations in the stability 
derivatives with respect to u used in the three -degree -of -freedom dynamic analysis. 

To illustrate the significance of this feature, stability analyses were made for a single 
(forward) c.g. location with the dynamic -pressure contribution to the derivatives X(i, 

Zq, and omitted. The aeroelastically corrected Mach number contributions were 

included, however. The resulting characteristic stability roots are compared in figure 17 
with those from the more complete analyses under the assumption of a locally uniform 
atmosphere. The omission of the q contribution had little effect on the short -period 
roots except for those cases in which coupling between the short -period and phugoid roots 
resulted in a coupled oscillatory mode. The omission had a strong destabilizing effect on 
the phugoid mode, however, particularly for the subsonic conditions, and resulted in an 
aperiodic divergence for all but the lowest dynamic pressures. The resulting effect on 
the static control parameters is shown in figure 18 by the curves faired with solid lines. 
As would be expected from the characteristic roots, omission of the q contribution 
had little effect on the maneuvering control parameter 5p|np but resulted in strongly 
unstable values of the speed control parameter 6p^u. 

Alternative Formulations for Normal-Acceleration Contribution 

In the present formulation, the effects of aeroelastic distortion due to normal accel- 
eration are embodied implicitly in the stability derivatives as indicated by the terms in 
equations (18) involving partial derivatives with respect to n. One alternative formula- 
tion which has been used in some segments of the industry results from omitting these 
terms from the stability derivatives defined by equations (18) and retaining derivatives 
with respect to n as explicit coupling terms in the three force and moment equations 
of motion. The set of equations of motion must then include an additional equation such 
as equation (34) to define the role of these coupling terms in the dynamic solution. The 
resulting solution is identical to that from the present formulation, but approximate 
expressions for static stability and control parameters would appear more complicated 
than those presented herein. 
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Another alternative formulation is that denoted in reference 2 as "Formulation I" 
and in reference 3 as the "Direct Formulation." In these formulations, the normal- 
acceleration contributions are broken down into inertial derivatives with respect to 
variables representing linear acceleration w, centrifugal acceleration q, and gravity 
orientation Q. These inertial derivatives appear in the equations of motion along with 
aerodynamic derivatives representative of the "zero mass" elastic airplane. The 
dynamic solution of these equations is equivalent in all important respects to that from 
the present formulation with the reservation that the contribution of dynamic-pressure 
perturbations to the forward-speed derivatives is not included in the dynamic formula- 
tions of references 2 and 3. Although the inertial derivatives appear in combination with 
aerodynamic derivatives in the equations of motion, it is not recognized in references 2 
and 3 that these combined derivatives can be considered as stability derivatives generally 
applicable to static stability, control, and performance analysis, as well as dynamic 
analysis. 

The final formulation to be discussed is that denoted in reference 2 as "Formula- 
tion n" and in reference 3 as the "Indirect Formulation," It is stated in those references 
that this formulation "represents what has been the standard method of accounting for 
inertial effects, or as it is frequently called, inertial relief." The derivation of this for- 
mulation makes use of the expression 


n=l+np=l + 



(61) 


Note that this expression implies that n is the load factor rather than the component of 
acceleration normal to a structural reference plane. Note also that this expression is 
valid only for constant dynamic pressure. 

To describe the concept of the indirect formulation in the notation of the present 

paper. 


V 9Ct 

= ^L, 1 2/ "8^ 

k=l ^ 

(62a) 

K 

V 9C.V, 

Cm = Cm,l + 2 Pp,kl^ 
k=l ^ 

(62b) 

k=l ^ 

(62c) 
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where Pj^ is a set of K physical variables having n as a member. The indirect for- 
mulation makes use of equation (61) to show that the following expressions are equivalent 
to equations (62): 


K-1 j 

y p , 

/ aCL \ 

®Pk 1 

L Pp,k 

lr-1 \ 

1 ._JL_££Li 

Jv— X y 

^L,l 

K-1 

/ 


Cm - Cm, 1 + ^ Pp,k 
k=l 


8C™ apj^ 


% aCL 

Cl,i -“aH- 


V 


K-1 


Cd = ^D,1+ ^ Pp,l 


k=l 


( 9Cp 

9Cp an apj^ 


(63a) 


(63b) 


(63c) 


where n is not included in the set of K-1 physical variables Pj^. The expressions 
in parentheses in equations (63) can then be transformed into stability derivatives with 
respect to the motion variables by a procedure analogous to that used in the present 
formulation. 

The indirect formulation was applied to the forward c.g., uniform -atmosphere cases 
of the present analysis. The resulting characteristic roots are compared with those from 
the present formulation in figure 19, and the static control parameters obtained from the 
transfer function solutions are shown in figure 18. Although the contributions of dynamic - 
pressure perturbation were not included in the indirect formulation of references 2 and 3, 
it is clear from equations (63) that q could be either included in or excluded from the 
set of physical variables. Both conditions are represented by the results given in fig- 
ures 18 and 19. The effects of omitting the q contribution in the indirect formulation 
are similar to the effects already discussed relative to the present formulation. 

The indirect formulation predicts a significantly greater speed stability than the 
present formulation as evidenced by larger values of both the speed control parameter 
6p/u ^d the phugoid frequency. This discrepancy of speed-related characteristics 
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should be anticipated from the invalidity of equation (61) imder conditions of varying 
dynamic pressure. A more subtle inadequacy of the indirect formulation is revealed by 
observing that although the maneuver control parameter 6p/np is accurately predicted^ 
the magnitude of the short -period roots is consistently underpredicted. This discrepancy 
is apparently related to the inappropriate application of inertial relief, imder the indirect 
formulation, to all derivatives including those describing control effectiveness. 

CONCLUDING REMARKS 

A set of longitudinal stability derivatives has been defined, and the corresponding 
equations of motion formulated, applicable to the analysis of longitudinal static and 
dynamic stability, control, and performance of elastic airplanes. This development is 
subject to the assumptions of small perturbations from a steady reference flight condition, 
structural deflection constrained to a direction normal to a structural reference plane, 
and the quasi-steady aeroelastic assumption of structural deflections proportional to 
applied loads. The development described herein avoids any constraint on the forward- 
speed degree of freedom, and the restilting stability derivatives exhibit two important 
departures from past practice: (1) The aeroelastic contributions of dynamic -pressure 

perturbations are included, and (2) the aeroelastic contributions of normal acceleration 
appear primarily in the derivatives with respect to pitching velocity and angle -of -attack 
rate and, for an imaccelerated reference flight condition, do not influence the derivatives 
with respect to angle of attack. 

Approximate expressions for the static and maneuvering stability margin and con- 
trol parameters in terms of the stability derivatives defined herein have been shown by 
means of illustrative calculations to correlate well with the characteristic stability roots 
and response to control inputs calculated by solution of the dynamic equations of motion. 
Some results of illustrative analyses of an elastic airplane cast doubt on the general exis- 
tence of a neutral point or maneuver point where these points are defined as center -of- 
gravity locations corresponding to zero static margin or maneuver margin, respectively. 
Some of the results for an elastic airplane at high dynamic pressure indicate that for 
stable values of both the static margin and maneuver margin, as the maneuver margin 
becomes small, coupling between the short -period and the phugoid modes can result in a 
strongly unstable phugoid oscillation at a significantly increased frequency. Further 
increases in phugoid frequency result from consideration of the variation of atmospheric 
properties with altitude. 

The present formulation of the longitudinal stability and control problem for elastic 
airplanes is shown to be more complete than several alternative formulations, and one 
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commonly used alternative is shown to yield significantly inaccurate predictions of not 
only the speed stability and phugoid characteristics, but also the short-period frequency 
and damping. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 21, 1972. 
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APPENDIX 

A FORTRAN PROGRAM FOR CALCULATING LONGITUDINAL STABILITY 
DERIVATIVES OF ELASTIC OR RIGID AIRPLANES 

The FORTRAN program described herein calculates the longitudinal stability deriv- 
atives of an elastic airplane using the procedure outlined in the section of this paper enti- 
tled "Calculation of Stability Derivatives." The aerodynamic and structural influence - 
coefficient matrices and the airplane jig shape are not calculated in the program but must 
be supplied as inputs. The program is written in CDC FORTRAN IV, version 2.3 to run 
on Control Data 6000 series computers with the SCOPE 3.0 operating system and library 
tape. The statements in the program most likely to require change to run on other oper- 
ating systems are those calling the library subroutine MATRIX, which is used in this 
program only for matrix multiplication and inversion. The program is dimensioned to 
accommodate up to 120 panels on one side of the aircraft plane of symmetry and requires 
a field length of approximately 140 OOOg words. Input and output data are in U.S. Cus- 
tomary Units. 


Input Data 

The force sign convention used in the program is such that positive forces act in 
the negative Z -direction. The elements of both the aerodynamic and the structural 
matrices, therefore, have signs opposite to those of the and [s] matrices as 

defined in the body of this paper. Units of the matrix elements are feet^ for the 
matrix and 1/pound for the [s] matrix. 

Tape inputs . - TAPE 1 must contain the elements of three aerodynamic influence - 
coefficient matrices: The first corresponds to the nominal Mach number AM, the sec- 
ond to a slightly higher Mach number AM + DM, and the third to a slightly lower Mach 
number AM - DM. A value of DM of about 2 to 4 percent of AM is suggested. Each 
matrix has dimensions M x M, where M is the number of panels on one side of the air- 
craft plane of symmetry. The matrix elements are arranged in the following order: 

A(l, 1), A(l,2), . . ., A(1,M), A(2, 1), . . ., A(2,M), . . ., A(M,M), where the first index is 
the aerodynamic control point (slope point) index and the second is the load point index. 

No end-of-file code is to be used between successive matrices. 

TAPE 2 must contain the elements of two structural slope influence -coefficient 
matrices: The first relates slopes at slope points to unit loads at load points, and the 

second relates slopes at load points to imit loads at load points. Each matrix has dimen- 
sions M X M, and the matrix elements are arranged in order: S(l, 1), S(l,2), . . ., 

S(2,M), . . ., S(M,M), where the first index is the slope point index and the second is the 
load point index. No end-of-file code is to be used between the two matrices. 
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APPENDIX — Continued 


Card inputs. - The program is organized to execute multiple cases by use of nested 
DO loops. The outer loop provides for various center -of -gravity locations, and succes- 
sive inner loops provide for various airplane weights at each center of gravity, various 
dynamic pressures at each weight, and various load factors at each dynamic pressure. 
Level 1 input cards apply to all cases and are read once. Level 2 cards must be supplied 
for each center of gravity, level 3 for each weight and center of gravity, and level 4 for 
each dynamic pressure, weight, and center of gravity. 

Level 1 

Card 1, Format(3Il,2I3,5F10.0) 

IFLEX Set 1 for elastic case, requires TAPE 1 and TAPE 2. Set 0 for rigid case, 
requires TAPE 1 only. 

LIST Set 1 for output listing of elastic air loads and flexible slope increments 

due to jig shape, alpha, delta, and normal acceleration and total slopes 
for reference flight condition. Set 0 to delete above listings from output. 

LFT Set 1 for banked -turn reference flight condition. Set 0 for wings -level 

reference flight condition. 

M number of panels on one side of aircraft plane of symmetry 

NC number of center -of -gravity locations 

AM Mach number of reference flight condition 

DM Mach number increment used in preparing TAPE 1 

SREF reference wing area, feet^ 

CBAR reference chord, feet 

CAF increment in axial -force coefficient (skin friction, wave drag due to volume, 

etc.) to be added to that calculated from pressures normal to mean cam- 
ber surface. 

Each of the following arrays is required in format (8F10.0), I ranges from 1 to M: 
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APPENDIX — Continued 

TSR(I,1) mean-camber -surface slope at control point I due to jig shape, radians 

TSR(I,3) increment in mean-camber -surface slope at control point I due to unit 6, 

radians 

TLRJ(I) mean-camber -surface slope at load point I due to jig shape, radians 

TLRD(I) increment in mean -camber -surface slope at load point I due to unit 6, 

{'^6,f}, radians 

XSO(I) x-coordinate of slope point I from nominal origin, X positive forward, feet 

XLO(I) x-coordinate of load point I, feet 

Level 2 

Card 1, Format(I3, FIO.O) 

NW number of airplane weights 

XCG x-coordinate of center of gravity from nominal origin, X positive forward, 

feet 

Level 3 

Card 1, Format(I3,4F10.0) 

NQ number of dynamic pressures 

W airplane weight, pounds 

YI moment of inertia in pitch, ly, slug-feet^ 

Y1 (I 2 - IxY^Y value of 1.0 is appropriate if the actual value is unknown.) 

Y2 %z/^Y value of 0, is appropriate if the actual value is unknown.) 

The following array is required only if IFLEX has been set to 1, format (6F12.0), I ranges 
from 1 to M: 

FI(I) weight assigned to panel I, assumed to be located at load point I, poimds 
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Level 4 

Card 1, Format(I3,6F10.0) 

NN number of load-factor values 

Q dynamic pressure in reference flight condition, pounds/foot^ 

V true airplane velocity in reference flight condition, feet/second 

G gravity acceleration in reference flight condition, feet/second^ 

RHO air density in reference flight condition, Pj, slugs/foot^ 

RHOH (^jpV 1/foot 

Card(s) 2, Format(8F10.0) 

FL(I) load factor in reference flight condition, I ranges from 1 to NN (Note that 

FL = G/32.174 for a straight, level reference condition.) 

Output Data 

Program output is listed in two parts. The first part is identified by the values of 
Mach number, center -of -gravity location, weight, and dynamic pressure read in as input. 
This part presents a tabulation over all panel load points of the aerodynamic normal force 
divided by dynamic pressure and the flexible increment in mean -camber -surface slope 
produced by the jig shape and by unit values of angle of attack, control deflection, and 
normal acceleration. 

The second part is repeated for each input value of load factor and presents data in 
five groups identified as trim values, derivative contributions, stability derivatives, static 
parameters, and total trim condition slopes. The trim values define the reference flight 
condition and include 1 , , 1j ^m, 1j ^ 1 , ( p , nj, and qj^cy’2Vi along 

with the number of iterations used to solve the nonlinear reference equations. The deriv- 
ative contributions include the increments in Ci«j and due to jig shape and the par- 

tial derivatives of Ca, Cn, and Cm with respect to a, 6, qc/2V, n, q, M, and q. 
The stability derivatives are presented in both coefficient form (Ca, Cn, and Cm) and 
dimensional form (X, Z, and M) as derivatives of these quantities with respect to u, u. 
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a, oi, 0, q (or qc/2V in the case of the coefficient derivatives), q, 6, and h. 
Unsteady aerodynamic effects are not considered in the program, and therefore the val- 
ues of the derivatives with respect to u, a, and q reflect only the aeroelastic effects 
arising from inertial loading. The static parameters include the stability derivative 
ratio values of static margin, maneuver margin, 6p/u, and dp^np 

calculated from equations (59) and (60) of this paper. The mean-camber -surface slopes 
in the reference flight condition are tabulated over all panel slope points and load points. 

The imits of all output quantities not defined as input data are as defined in the main 
text of this paper except that normal force per unit dynamic pressure is in feet^ and alti- 
tude is in feet. 


Sample Cases 

Listings of input data cards and program output for two sample cases are given on 
the following pages. Both cases correspond to the airplane used for the illustrative 
analysis of this paper flying at its supersonic design condition. Sample case 1 is for 
the rigid airplane, and sample case 2 is for the flexible airplane. The nominal origin 
of x-coordinates is at 0.33c. The airplane is represented by 110 panels numbered con- 
secutively from leading edge to trailing edge in rows progressing from root to tip. The 
supersonic paneling arrangement includes 11 chordwise rows with 10 panels per row. 
Data tabulated over all panels are listed in the order of panel numbers. 

Program compilation and execution required 7 seconds of central processor time 
for sample case 1 and 48 seconds for sample case 2. 
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INPUT DATA CAROS - SAPPLE CASE 1 


COLUMN NUPBER 

COOCOOC001inilllll22222222223333333333AA^4AAAA^45555555555666666666677777777778 

123456709O123A56709O123456789O123456789O123A5678SOX2345678<3C123456789O123A56789O 


0101 10 12.7 .1 9932.37 112.8618 .0C533 


-.0020858 

.0217330 

.0649896 

.0721554 

.0677656 

.0584697 

.0391723 

•0218184 

-.0198222 

-.0640468 

-.0173102 

.0054653 

.0166327 

.0223163 

.0250962 

.0253127 

.0236214 

.0206650 

.0164318 

•0115639 

-.0407757 

-.0091252 

-.0024986 

.0047684 

.0093366 

,0117246 

•0126243 

•0128901 

•011(661 

.0102973 

-.0523926 

-.0241893 

-.0143258 

-.0061958 

-.0006557 

.0026021 

.0053965 

.0064772 

.0073860 

.0080724 

-.0646976 

-.0367379 

-.0204279 

-.0144045 

-.0089829 

-.0049710 

-.0015805 

.0008318 

.0026083 

•0040618 

-.0768725 

-.0474848 

-.0314948 

-.0236675 

-.0175863 

-.0126706 

-.0091525 

-.0057909 

-.0031599 

-.0009936 

-.0051170 

-.0639445 

-.0418934 

-.0302575 

-.0242494 

--.02C8992 

-.0176371 

-.0140687 

-.0112759 

-.0093210 

-.0979577 

-.C725046 

-.0571262 

-.0479308 

-.0421487 

-.0357561 

-.0312538 

-.0292489 

-.0267000 

-.0228283 

-.0666034 

-.0616627 

-.0569024 

-.0507564 

-.0464143 

-.0432134 

-.0397698 

-.0373008 

-.0355413 

-.0340052 

-.0292534 

-.0403315 

-.0423783 

-.0430595 

-.0435389 

-.0438973 

-.0442617 

-.0445019 

-.0446310 

-.0446310 

.0063823 

-.0255481 

-.0370163 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 




0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0, 

0. 

0. 

0. 

0. 

0. 

0. 

c. 

c. 

0. 

0. 

0. 

C. 

0. 

0. 

0. 

.0532 

.459 

0. 

0. 

0. 

0. 

0. 

c. 

C. 

0. 

.3725 

.734 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

.4627 

.826 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

c. 

.049 

1 .0 

0. 

0. 

C. 

0. 

0. 

0. 

0. 

,2738 

.758 

.758 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

c. 

c. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

c. 

0. 




-.0436609 .1321082 .0554699 .C755199 .C738231 .C657603 .0510243 .0254630 

.0046529 -.0542514 -.0404669 -.0031124 .0133642 .0208076 ,0247195 .0260645 

.0248230 .0224995 .0186951 .01 37779 C6 15C 84 -.0174436 -.0048858 .0023824 

.0079625 .ClllOOO .0127531 .0132559 .C126467 .0110400 -.0688493 -.0334654 

-.0183140 -.CC89998 -.0021476 .0012953 .C05C137 .0060814 .0071681 .0078364 

-.0666560 -.0458547 -.0260490 -.0171017 -.0105445 -.0065365 -.0026051 .0000706 

.0020907 .0035243 -.0961399 -.0585720 -.0374848 -.0260449 -.0201491 -.0140063 

-.0105627 -.0069518 -.0039608 -.0018017 -.C883117 -.0722898 -.0505014 -.0335301 

-.0255958 -.0219175 -.0189672 -.0153701 -.0118629 -.0099270 -.0844228 -.0872726 

-.0613408 -.0487982 -.0435496 -.0370947 -.C30E534 -.0283116 -.0269155 -.0230386 

-.0637906 -.0650822 -.0602126 -.0542942 -.0460599 -.0454076 -.0415073 -.0383499 

-.0363429 -.0343914 -.0176354 -.0369715 -.C42C562 -.0428395 -.0434637 -.0437361 

-.0441703 -.0444051 -.0445780 -.0446094 .0253262 -.0117544 -.0343573 -.0437926 



04366C9 -. 

C436609 -.0436609 -.0436609 -.0436609 -.0436609 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0, 

0. 

0. 

C. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

C. 

0. 

0, 

0. 

0. 

c. 

0. 

0. 

.0532 

.459 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

.3725 

.734 

0. 

c. 

0. 

0. 

0 . 

0. 

0. 

0. 

G. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

.4627 

.826 

c. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

.849 

1 .0 

c. 

0. 

0. 

0. 

0. 

C. 

C. 

.2738 

.758 

.758 

0. 

c. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

c. 

0. 

0. 

0. 

c. 

0. 

0. 

0* 

c. 

0. 








dea 




> Edesff 




6,f 


Level 1 
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106.716 

84.264 

61.810 

39.356 

16.901 

-5.553 

-28.007 

-50.461 




-72.515 

-55.372 

59.089 

44.755 

30.422 

16.089 

1.756 

-12.578 




-26.511 

-41.244 

-55.577 

-69.914 

40.796 

28.307 

15.977 

3.568 




-8.841 

-21.250 

-33.660 

-46.C69 

-56.476 

-70.808 

22.414 

11.835 




1.256 

-5.323 

-19.901 

-30.480 

-41.059 

-51.638 

-62.217 

-72.796 




3.947 

-4.892 

-13.730 

-22.569 

-31.407 

-40.246 

-49.084 

-57.922 




-66.761 

-75.599 

-14.613 

-21.809 

-29.C04 

-36.200 

-43.395 

-50.591 

, 



-57.766 

-64.982 

-72.177 

-79.373 

-33.227 

-38.835 

-44.443 

-50.052 

y xa 



-55.660 

-61.269 

-66.877 

-72.486 

-78.094 

-83.703 

-51.854 

-55.890 




-55.526 

-63.962 

-67.998 

-72.C35 

-76.C71 

-80.IC7 

-84.143 

-36.175 




-66.429 

-65.252 

-72.C75 

-74.898 

-77.722 

-80.545 

-83.368 

-86.191 




-39.C34 

-51.837 

-76.963 

-78.944 

-8C.526 

-82.9C7 

-84.889 

-86.871 




-88.652 

-5C.834 

-92.816 

-54. 757 

-87.514 

-88.672 

-89.831 

-9C.909 




-92.148 

-93.306 

-54.465 

-95.623 

-56.782 

-57.94C 





y Level 1 

121.237 

98.029 

74.820 

51.612 

28.403 

5.195 

-18.014 

-41.223 




-64.421 

-87.640 

65.759 

51.403 

37.C47 

22.651 

8.335 

-6.C21 




-20.378 

-34.734 

-45.C90 

-63.448 

46.620 

34.138 

21.755 

9.322 




-3.110 

-15.543 

-27.975 

-40.408 

-52.841 

-65.273 

27.446 

16.842 




6.238 

-4.366 

-14.971 

-25.575 

-36.179 

-46.783 

-57.388 

-67.992 




6.230 

-.635 

-5.500 

-18.366 

-27.231 

-36.096 

-44.961 

-53.826 




-62.651 

-71.557 

-11.017 

-18.242 

-25.468 

-32.693 

-39.918 

-47.143 




-54.368 

-61.594 

-68.019 

-76.C44 

-3C.251 

-35.896 

-41.541 

-47.187 




-52.822 

-58.477 

-64.122 

-65.767 

-75.413 

-81.050 

-49.410 

-53.497 




-57.564 

-61.671 

-65.759 

-69.846 

-73.933 

-78.020 

-32.107 

-86.195 




-64.884 

-67.729 

-7C.573 

-73.418 

-76.262 

-79.1C7 

-81.952 

-84.796 




-87.641 

-50.486 

-75.690 

-77.701 

-75.711 

-81.722 

-83.732 

-85.743 




-87.753 

-89.764 

-91.775 

-93.765 

-06.35C 

-87.557 

-88.763 

-85.970 




-51.176 

-52.383 

-93.599 

-94.796 

-96.002 

-57.209 


J 




1-9.823 










Level 2 

1^72506. 

4256C516. 

1. 

0. 







Level 3 

1 59 6 . 6 1 5 

2613.62 

31 .573 

.CCC17A65 

-.CCCC475 

C. 




1 

^ Level h 

993753 
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PROGRAM OUTPUT - SAMPLE CASE 1 


AEROELASTIC DERIVATIVES, M= 2.700, XCG= -9. 


ELASTIC AIRLOADS 

PER Q DUE TC 

JIG SHAPE 

4.405167 

5.086942 

6.407287 

5.500033 

-8.101925 

-5.825629 

2.066525 

3.56O051 

.371886 

-.545671 

-1.497147 

-2.352072 

2.124071 

1.518737 

1.015405 

.307909 

1.846014 

2,000352 

1.950671 

1.554878 

-3.577497 

.861491 

1.613613 

1.468820 

.749346 

.489162 

-4.567554 

.003242 

1,157658 

.988819 

.800750 

.650775 

.794494 

.7C7317 

.726670 

.774167 

-.3C 6045 

.263C62 

.378186 

.370769 

-3.965022 

-2.247014 

-1.251506 

-.588971 

.162971 

.085693 

-1.451682 

-2.063045 

-.522849 

-.289209 

-.119487 

-.023573 

-1.165501 

-1.117301 

-1.034143 

-.942978 

ELASTIC AIRLOADS 

PER Q DUE TC 

ALPHA 

149.810868 

99.252230 

99.386244 

116.837692 

104.610133 

52.649335 

170.638300 

92.936603 

75.136618 

74.665779 

72.390064 

66.481118 

71.263026 

70.013380 

68.565508 

68.104686 

82.799433 

72.631784 

67.249463 

65.149742 

182.627076 

102.495503 

81.837039 

70.107545 

54.687695 

53.506766 

177.915649 

100.487495 

52.180692 

50.892719 

49.937336 

48.327395 

55.441732 

51.305983 

47.478609 

44.581055 

65.643622 

52.061183 

45.853602 

45.481612 

109.312972 

78,634141 

60.797135 

49.212877 

30.735940 

30.495578 

58.849529 

64.530916 

33.034518 

28.6569C7 

25-071517 

22.740085 

33.654015 

32.852887 

31.404300 

29.787274 


823, W= 472500, QBAR = 

596.62 


4.137748 

1.658043 

-1.973766 

-4.281655 

3.310008 

2.858762 

2.303882 

1.333546 

-.653648 

3.177942 

2.440543 

2.408191 

-.344136 

-.949054 

-2.276323 

1.950561 

1.197829 

.774787 

.332770 

-.040082 

1.537863 

1.432235 

1.180391 

.962500 

1.076597 

1.013636 

1.113711 

1.216811 

-5.144581 

-1.C41628 

.626273 

.865549 

.737781 

.6008C2 

-5.457209 

-1.526993 

.320750 

.242784 

.295230 

.434817 

-.15C487 

. 135972 

.244314 

.231376 

-1.750670 

-1.370496 

-1.039768 

-.781034 

.C55301 

-.623199 

-.926234 

-1.128723 

-.849150 

-.752640 



126.397336 

134.199253 

137.049141 

134.909163 

78.378010 

75.325785 

74.677153 

75.125170 

180.308611 

100.213469 

81.761434 

74.373917 

66.901651 

64.536145 

183.797293 

102.669730 

63.293960 

61.576594 

6C. 865743 

59.815470 

63.472436 

59.577503 

58.183180 

56.485101 

77.388528 

66.896804 

59.962141 

55.245622 

169.214616 

95.769864 

71.727225 

61.029331 

43.062685 

43.184773 

155.538850 

08.224507 

44.911300 

43.635639 

41.604817 

39.146471 

41.186243 

34. 759944 

32.125564 

31.094826 

57.373415 

49. 801898 

43.254497 

38.002822 

26.594636 

28.5840 89 

30.718J88 

32.887397 

28.086723 

26.314033 




ELASTIC AIRLOADS PER Q DUE TO DELTA 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

O.ODOOOO 

0.000000 

.410456 

3.796512 

5.776272 

0.000000 

0-000000 

0.000000 

0.000000 

0.000030 

0.000000 

O.OOOCOO 

.036935 

.703706 

8.182117 

o.ocoooo 

0.000000 

0.000000 

0.000000 

0.000000 

O.OOOCOO 

0.000000 

.008863 

5,904510 

53.047485 

0.000000 

0,000000 

0.000000 

o.oocooo 

0.000000 

0.000000 

0.000000 

.001893 

32.852237 

50.013243 

O.OOOCOO 

0,000000 

0.000000 

O.OOOOQC 

0.000000 

0.000000 

0.000000 

.000347 

4.269627 

18.697504 

o.oocooo 

0.000000 

o.ocoooo 

0.000000 

0.000000 

O.OCOOOO 

0.000000 

.000053 

31.418286 

48.231617 

0,000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

O.OPOOOO 

.000005 

49. 519799 

52.151664 

0.000000 

0.000000 

0.000000 

o.acoooo 

0.000000 

0.000000 

0.000000 

12.447534 

30.458912 

30.913052 

0.000000 

0.000000 

0.000000 

0.000000 

0,000000 

O.OODDOO 

0,000030 

0.000000 

.596008 

2.451511 

0.000000 

0.000000 

0.000000 

O.OOOOOn 

0.000000 

o.orooco 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0-000000 

0.000000 

0.000000 

0.000000 

O.OOOCOO 

o.cooooo 

0.000000 

o.oocooo 

0.000000 
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APPENDIX - Continued 


LOAD FACTOR « *99375 


TRIM VALUES ( 2 ITERATIONS) 


CA * 

.004759 

ALPHA - 

•048506 

PHI 

0.000000 

CN » 

.079145 

DELTA - 

•000002 

AN* G - 

•992584 

CM » 

0.000000 

THETA » 

.048506 

QC/2V - 

.000000 


DERIVATIVE CONTRIBUTIONS 



JIG SHAPE 

ALPHA 

DELTA 

QC/2V 

AN* G 

QDOT 

MACH 

OYN PRESS 


CA 


-.0301833 

.0031426 

.0038041 

0.0000000 

0.0000000 

•0006716 

0.00000000 

CN 

.003698 

1.555408 

•088980 

•672663 

Q. 000000 

0.000000 

-.018491 

0.0000000 


CM 

•009699 

-.199950 

-.049690 

-.539630 

0.000000 

0.000000 

•003161 

0.0000000 



STABILITY 

DERIVATIVES 









U 

UDQT 

ALPHA 

ALPHA DOT 

THETA 

Q 

QDOT 

DELTA 

ALTITUDE 

CA 

.0018134 

-0.0000000 

-.0301833 

-0.0000000 

-0.0000000 

•0038041 

0.0000000 

•0031426 

0. 

CN 

-.049926 

-0.000000 

1.555406 

-0.000000 

-0.000000 

•672663 

0.000000 

•088980 

0. 

CM 

•008535 

-0.000000 

-.199950 

-0.000000 

-0.000000 

-.539630 

0.000000 

-.049690 

0. 

X 

-.001749 

-0.000000 

.004660 

0.000000 

-0.000000 

-.000013 

-0.000000 

-.000485 

3.48992E-08 

Z 

-.016729 

-0.000000 

-.240115 

0.000000 

-0. 000000 

-.002242 

-0.000000 

-.013736 

5.80350E-07 

M 

•134113 

0.000000 

-3.141995 

0.000000 

0.000000 

-.183071 

0.000000 

-.780820 

0. 


STATIC PARAMETERS 


CH ALPHA/CN ALPHA 
STATIC MARGIN 
MANEUVER MARGIN 
DELTA/U 
OELTA/ANf G 


••128552 

-.141924 

-.130062 

•587300 

-.271118 


TOTAL TRIM CONDITION SLOPES AT SLOPE POINTS 


-.0020858 

.0217330 

.0649896 

.0721554 

.0677656 

•0584697 

.0391723 

•0218184 

-.0198222 

-.0640468 

-.0173182 

.0054653 

•0166327 

•0223163 

•0250962 

•0253127 

.0236214 

.0206850 

.0164318 

•0115639 

-.0407757 

-.0091252 

-.0024986 

•0047684 

.0093366 

.0117246 

•0126243 

.0128901 

.0118662 

.0102982 

-.0523926 

-.0241893 

-.0143258 

-.0061958 

-.0008557 

•0028021 

.0053965 

.0064772 

.0073867 

•0080739 

-.0646976 

-.0367379 

-.0204279 

-.0144045 

-.0089829 

-.0049710 

-.0015805 

•0008318 

•0026083 

•0040618 

-.0768725 

-.0474848 

-.0314948 

-.0236675 

-.0175863 

-.0126706 

-.0091525 

-.0057909 

-.0031590 

-.0009919 

-.0851170 

-.0639445 

-.0418934 

-.0302575 

-.0242494 

-.0208992 

-•0176371 

-.0140687 

-.0112742 

-.0093190 

-.0979577 

-.0725046 

-.0571262 

-.0479388 

-.0421487 

-.0357561 

-.0312538 

-.0292484 

-.0266985 

-.0228268 

-.0668034 

-.0616627 

-.0569024 

-.0507564 

-.0464143 

-.0432134 

-.0397698 

-.0373008 

-.0355413 

-.0340052 

-.0292534 

-.0403315 

-.0423783 

-.0430595 

-.0435389 

-.0438973 

-.0442617 

-.C445019 

-.0446310 

-.0446310 

•0063823 

-.0255481 

-.0370163 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 



TOTAL 

TRIM CCNOITICN SLOPES 

AT LOAD POINTS 




-.0436609 

•1331082 

.0554699 

.0755199 

.0738231 

.0657603 

•0510243 

•0254630 

•0046529 

-.0542514 

-.0404669 

-.0031124 

.0133642 

• 020 8076 

.0247195 

•0260645 

.0246230 

•0224995 

.0186951 

.0137779 

-.0615084 

-.0174436 

-.0048858 

.0023824 

•0079625 

•0111000 

•0127531 

.0132559 

.0126468 

.0110417 

-.0688493 

-.0334654 

-.0183140 

-.C089998 

-.0021476 

.0012953 

.0050137 

•0060814 

•0071688 

.0078379 

-.0866568 

-.0458547 

-.0260490 

-.0171017 

-.0105445 

-.0065365 

-.0026051 

•0000706 

.0020 907 

.0035243 

-.0961399 

-.0585720 

-.0374848 

-.0260449 

-.0201491 

-.0140063 

-.0105827 

-.0069518 

-.0039599 

-.0018000 

-.0883117 

-.0722898 

-.0505014 

-.0335301 

-.0255958 

-.0219175 

-.0189672 

-.0153701 

-.0118612 

-.0099250 

-.0844228 

-.0872726 

-.0613408 

-.0487982 

-.0435496 

-.0370947 

-.0308534 

-.0283111 

-.0269140 

-.0230371 

-.0637906 

-.0650822 

-.0602126 

-.0542942 

-.0480599 

-.0454076 

-.0415073 

-.0383499 

-.0363429 

-.0343914 

-.0176354 

-.0369715 

-.0420562 

-.0428395 

-.0434637 

-.0437361 

-.0441703 

-.0444051 

-.0445780 

-.0446094 

.0253262 

-.0117544 

-.0343573 

-.0437926 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 

-.0436609 
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APPENDIX - Continued 


INPUT DATA CAPOS - SAMPLE CASE 2 


COLUMN NUMBER 

0OO000C0OllimillX2222222222333333 32 33444^44444455 55 555555666666666677777777778 
123^5678SC123A5678901234567890123A56789C123A5678S0123A56785C1234567890123456789C 


IIOIIO 12.7 

.1 

9932.37 

112.8618 

.00533 



> 


.0015623 

.C246630 

.06678C2 

.C729176 

.0678190 

.0584583 

.0388353 

.0199932 



-.0246213 

-.C674348 

-.0162212 

.0062563 

.0170591 

.0224678 

.0250616 

.0253794 



.0235869 

.C2C7695 

.0244825 

.CC93577 

-.C410C27 

-.0091921 

-.0027604 

.0C43669 



.009C743 

.C118738 

.0133369 

.0180953 

.0172433 

.0153694 

-.0533342 

-.0245784 



-.0143609 

-.0055302 

.0006991 

.0057368 

.0096975 

.0171323 

.0X95188 

.0200882 



-.0657932 

-.C37C484 

-.C191335 

-.0111876 

-.0034572 

.0036071 

.0101099 

.0174297 



.C2065C8 

.C224772 

-.0738053 

-.0426467 

-.0239943 

-.0141310 

-.0049521 

.CC37C86 


>■ G . . 

.0112772 

.C203119 

.0303789 

.0318917 

-.0738209 

-.0499947 

-.0247094 

-.0101866 


Jig 

.0903692 

.CC70773 

.C143163 

.0223143 

.0304250 

.0314286 

-.0737579 

-.0450515 



-.0258222 

-.0121348 

-.CC13292 

.0065629 

.0139206 

.0166269 

.0193716 

.0224362 



-.02501C6 

-.C17540C 

-.0004620 

.C0059C6 

.0077007 

.0125698 

.0238697 

.0374343 



.C3939E7 

.C347754 

.C287487 

.C1928C4 

.0206758 

. C223097 

-C246000 

.02656J4 



.0317445 

.C4C9611 

.0414466 

.0389605 

.0809498 

.0509308 

.0415179 

.C380638 



.0404353 

.0439335 

.C406 843 

.0415755 

.04381C6 

.0459071 





C. 

0. 

0. 

0. 

C. 

0. 

0. 

0. 



0. 

0. 

0. 

C. 

c. 

0. 

0. 

0. 



0. 

0. 

0. 

c. 

c . 

0 . 

0. 

0. 



0. 

0. 

0. 

0. 

.0532 

.459 

0. 

0. 



c. 

0. 

0. 

0. 

0. 

0. 

.3725 

.734 



c. 

0. 

0. 

0. 

C. 

C. 

0. 

0. 



0. 

c. 

c. 

0. 

0. 

0. 

0. 

0, 



0. 

0 . 

-4627 

.626 

c. 

c. 

0. 

0. 


6 

c. 

c. 

0. 

0. 

.849 

1 .0 

0. 

0. 



c. 

0. 

0. 

0. 

0. 

. 27 38 

.758 

.758 



c. 

0. 

0. 

0. 

C. 

0 . 

0. 

0. 



0. 

0 , 

0. 

0. 

C. 

0 . 

0. 

0. 



c. 

0. 

0. 

0. 

C. 

0. 

0. 

0. 



0. 

0. 

0. 

0. 

0. 

c. 





-.0395242 

. 1368794 

.C57968C 

.C768654 

-0741514 

.0657412 

.0509240 

.0247068 



.0010444 

-.C584475 

-.0392262 

-.0021710 

.0135810 

.0210751 

.0247217 

.0261119 



.0248561 

. C224076 

.C179596 

.0114259 

-.0619094 

-.0175373 

-.0049803 

.0020224 



.C076296 

.C111179 

.0125381 

.0168279 

.0170874 

.0162129 

-.0700422 

-.0340887 



-.C134C99 

-.CC66985 

-.CC11399 

.0035465 

,0079648 

.0144903 

.0188801 

.0199499 



-.0878‘:22 

-.0466441 

-.0256543 

-.0148723 

-.0062754 

.0005924 

.0073152 

.014807? 



.C1952C4 

.0217157 

-.C94C483 

-.0548035 

-.0314223 

-.01741 16 

-.0091469 

.OOC5C58 


L c- 

.00754EC 

.0159857 

.C2 70648 

.0319256 

-.0782119 

--C5988C1 

-.0350143 

-.0153127 


^ jig,: 

-.0035691 

.CC45208 

.C1C6698 

.0189262 

.0267759 

.0315315 

-.0615457 

-.0620748 



-.0324779 

-.0156838 

-.0047259 

.CC366CC 

.0123830 

.0173652 

.0186121 

.0221058 



-.0227201 

-.0226C67 

-.C157763 

-.0046018 

.004578C 

.0098240 

.0172427 

.0332366 



.033987C 

.C369C10 

.C390C59 

.0210111 

.0191081 

.0207521 

-02289^1 

-0255374 



.0285304 

.C369951 

.C410641 

.0399069 

.0985235 

.0629536 

.0422116 

.0358328 



.0382653 

.C419499 

.C438443 

.0412608 

.0429089 

. C448589 





C. 

C, 

0. 

0. 

0. 

0. 

0. 

0- 



C . 

0. 

0. 

C. 

0. 

C. 

C. 

0. 



C. 

C. 

0. 

C. 

0. 

0. 

0. 

0. 



0. 

C . 

0. 

0. 

. 0532 

.459 

0. 

0. 



c. 

c. 

0. 

0. 

C. 

0 . 

. 3725 

.734 



0. 

c • 

0. 

0. 

C. 

0 . 

0. 

0. 



c. 

c . 

0. 

0. 

C. 

0. 

0. 

0. 


^ E r -C* 

c. 

c . 

.4627 

.826 

0 . 

c. 

0. 

0. 


0 

0. 

0. 

0. 

0. 

.849 

1 .0 

0. 

0. 



0. 

c. 

0 . 

0. 

0- 

. 2738 

.758 

.758 



0. 

c . 

c. 

c. 

0. 

C. 

0. 

0. 



0. 

c . 

0. 

0. 

C . 

C . 

C. 

0. 



0. 

c. 

0. 

c- 

0. 

0. 

0. 

0. 



c. 

0. 

0. 

0. 

c . 

0 . 



J 
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106.718 e^.264 61.010 39.356 16.901 -5.553 -28.007 -50.461 

-72.915 -95.372 59.C89 44.755 30.422 16.089 1.756 -12.578 

-26.911 -41.244 -55.577 -69.914 4C.796 28.387 15.977 3.568 

-8.041 -21.250 -33.660 -46.069 -50.478 -70.888 22.414 11.835 

1.256 -9.323 -19.901 -30.48C -41.C59 -51.638 -62.217 -72.796 

3.947 -4.892 -13.730 -22.569 -31.407 -40.246 -49.084 -57.922 

-66.761 -75.599 -14.613 -21.6C9 -29.004 -36.200 -43.395 -50.591 

-57.766 -64.982 -72.177 -79.373 -33.227 -38.835 -44.443 -5C.C52 

-55.660 -61.269 -66.677 -72.486 -78.C94 -83.703 -51.854 -55.890 

-59.926 -63.962 -67.998 -72.035 -76.C71 -80‘.107 -84.143 -88.179 

-66.^29 -69.252 -72.075 -74.898 -77.722 -80.545 -83.368 -86.191 

-35.014 -91.837 -76.963 -78.944 -80.926 -82.907 -84.889 -86.871 

-88.852 -9C.834 -92.816 -94.797 -87.514 -88.672 -89.831 -9C.989 

-92.148 -93.3C6 -94.465 -95.623 -96.782 -97.940 

121.237 58.029 74.020 51.612 28.4C3 5.195 -18.014 -41.223 

-64.431 -67.640 65.759 51.403 37.C47 22.691 8.335 -6.021 

-20.378 -34.734 -49.090 -63.448 46.620 34.188 21.755 9.322 

• -3.110 -15.543 -27.975 -4C.408 -52.841 -65.273 27.446 16.842 

6.238 -4.366 -14.971 -25.575 -36.179 -46.783 -57.388 -67.992 

8.220 -.635 -9.500 -18.366 -27.231 -36.096 -44.961 -53.826 

-62.651 -71.557 -11.017 -13.242 -25.468 -32.693 -39.918 -47.143 

-54.368 -61.594 -68.819 -76.C44 -3C.251 -35.896 -41.541 -47.187 

-52.632 -56.477 -64.122 -69.767 -75.413 -81.058 -49.410 -52.497 

-57.564 -61.671 -65.759 -69.846 -73.933 -78.020 -82.107 -86.195 

-64.664 -67.728 -7C.573 -73.418 -76.262 -79.1C7 -81.952 -64.796 

-37.641 -90.486 -75.690 -77.701 -79.711 -81.722 -83.732 -85.743 

-37.753 -89.764 -91.775 -93.785 -86.35C -87.557 -88.763 -89.970 

-91.176 -52.363 -53.539 -54.796 -96.C02 -97.209 

1-9.623 

14725CC. 4256C516. 1. C. 

12392.932 5062.200 9545.400 8512.676 7968.411 8564.713 

1D4C8.682 31051.640 5265.812 7099.260 -566.222 5861.568 

273,140 7257.529 4890.585 5056. 2°0 1568.800 106-^.500 

5627. C76 3729. CIO 382.190 61,068 2540. 879 -518. 033 

llici. 306 1631. C30 1301.480 1285 . 46C 1553.850 8894,700 

6C3.32C 2833.442 3190.6C6 C.CCC 2915.319 2215.237 

22C6.C5C 2011.224 942.250 -23.882 664.260 1825.644 

1556.755 4428.328 2486.367 2418.258 2C93.0CC 1652.776 

1{45£.C17 7C4.51C 246.882 65C.899 871.566 1036.160 

1772.157 1188.848 986.766 975.430 1C86.C54 280.930 

593.273 1196.724 157.690 902.372 1126.448 852.943 

I4f4.395 -72.070 365. 1C6 336,080 779.323 664.872 

737.529 0.000 772.506 771.619 810.895 361.850 

485.940 375.130 72.122 67.310 67.575 293.290 

156.540 361.920 229.960 496.010 117.050 148.060 

202.010 0.000 112,950 145.890 -11,100 47,608 

17.619 56.599 2.417 68.295 127.410 42.469 

C.OCO 42.469 42.469 C.OOO 64.388 0.000 

64.366 0.000 

1596.615 2613.82 31.973 .C0C17465 -.CCC0475 0. 

.99 375 ? 
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APPENDIX - Continued 


PROGRAM OUTPUT - SAMPLE CASE 2 

AEROELASnC DERIVATIVES# 2.700, XCG* •9.823, W* 472500, QBAR= 596.62 
ELASTIC AIRLCAOS PER Q DUE TO JIG SHAPE 


5.810704 

5.793022 

6.845561 

5.747565 

3.965016 

.918621 

-2.722166 

-5,006404 

-8.567500 

-5.897203 

3.406742 

4.049156 

3.461013 

2.809673 

2.043070 

.927896 

.007622 

-.832394 

-1.798315 

-2.459552 

.422105 

3.613570 

2.557279 

2.324106 

1.897848 

1.281391 

.625798 

.345853 

-.391292 

-.912634 

-1.523217 

2.295400 

1.9^9128 

1.977447 

1.878833 

1.560983 

1.240003 

1.062064 

.56X957 

.198910 

-3.131752 

1*115642 

1.8371C1 

1.641616 

1.710963 

1.665462 

1.469050 

1.379527 

1.223425 

1.027146 

-4.106022 

.329446 

1.418165 

1.337896 

1.461160 

1.625356 

1.671620 

1.686506 

1.684915 

1.422207 

-4.430626 

-.539993 

1.087457 

1.333185 

1.357799 

1,299695 

1.423311 

1.565972 

1.597636 

1.412448 

-4.307457 

-.781343 

.332300 

,864164 

1.052227 

1.052091 

1.052806 

.982056 

1.030886 

1.159961 

-2.719483 

-1.313296 

-.440397 

.131023 

.521116 

.801604 

1.052415 

1.246255 

1.146631 

.898492 

-.279069 

-.954255 

-.749794 

-.464117 

--201243 

-.001218 

.265985 

.586256 

.740627 

.795548 

.990463 

.299702 

--014005 

-.211753 

-.259594 

-.229951 

-.229331 

-.183912 

-.123204 

-.056915 



ELASTIC AIRLOADS 

PER G DUE TQ 

1 ALPHA 





181.747655 

117.051915 

112-840567 

127.273047 

127.33C816 

102.307782 

92.614373 

81.607987 

53.749616 

23.664056 

204.221516 

106.640154 

85.393469 

75.638748 

67.436587 

53.893921 

48.112704 

44.413387 

34.523597 

29.612790 

211.412819 

112.278084 

86.477548 

71.116363 

59.323128 

46.249564 

41.394166 

26.882309 

24.5246C9 

24.953513 

208.707701 

111.283724 

83.712699 

64.915088 

51.245323 

41.223431 

34.472135 

19.153164 

17.554378 

20.617842 

198.649289 

1C6.C6E429 

76.167912 

55.925492 

43.213370 

33.141160 

27.311709 

18.638318 

15.975241 

15.799412 

177.6C8529 

94.736399 

65.662674 

48.966157 

37.180852 

28.768145 

22.431963 

15.686969 

10.85<547 

11.909068 

149.474245 

79.907197 

54.957269 

41,110378 

30.662899 

24.940629 

19.251603 

14.520721 

11.273539 

11.697674 

116.586145 

63,242091 

43.957051 

21.674819 

23.108955 

22.045206 

19.8C4528 

18.287747 

16.364706 

14.017691 

69.523591 

49.144325 

35.974603 

27.226434 

21.220555 

15.059433 

10.061113 

5.605167 

5.505273 

1.8967C1 

25-994001 

32.743521 

29.061866 

24-630426 

20.419141 

16.976281 

12.660722 

1.5830 80 

4.587620 

3.011404 

3.632053 

5.793406 

7.929713 

9.911498 

11.006C76 

1C. 865143 

11.142255 

10.682785 

9.925126 

9-040998 




ELASTIC AIRLOADS 

PER G DUE TO 

DELTA 





-.036512 

-.02C877 

-.015470 

-.000235 

.022378 

-1.405023 

-3.209775 

-4.672643 

-2.995317 

1.232063 

-.041329 

-.001117 

.059493 

-.021189 

-.204019 

-1.111716 

-2.062385 

-2.961556 

-4.583256 

2.124872 

.0*6355 

.068301 

.035234 

-.220115 

-.643060 

-1.3764C5 

-2.512645 

-5.943394 

-2.382481 

45.543862 

.073895 

.004073 

-.145799 

-.610768 

-1.232409 

-2.139722 

-3.4C0072 

-7.696791 

23.261931 

41.913720 

-.024656 

-.263321 

-.601959 

-1.199708 

-2.014875 

-3.534259 

-4.908C71 

-7.589171 

-3.671926 

10.658253 

-.630063 

-.788765 

-1.310599 

-2.180441 

-3,411694 

-4.663785 

-5.443196 

-5.773664 

21.49C796 

40.166445 

-2.472974 

-1.959083 

-2.237820 

-3.119007 

-4,557911 

-5.335321 

-4.450549 

-5.745179 

41.282109 

44.599598 

-5.380092 

-3.689807 

-3.314664 

-3.398184 

-4.297249 

-4.577133 

-4.967545 

7.598386 

25.752840 

26.242722 

-6.250481 

-4.5364T6 

-3.566472 

-3.057982 

-2.615758 

-2.610402 

-2.470375 

-2.312453 

-1.841749 

-.194833 

-4.313958 

-4.205368 

-2.543032 

-2.920926 

-2.416751 

-2.023792 

-1.632819 

-1.211640 

-1.1C5218 

-1-132324 

-1.637002 

-1.772055 

-1.837880 

-1.061547 

-1.815626 

-1.687032 

-1.609774 

-1.482136 

-1.329146 

-1.206799 




ELASTIC AIRLCACS 

PER C DUE TO 

AK, G 





-2.973795 

-1.58C326 

-1.C99015 

-.7598C6 

.IC7229 

2.296836 

2.920626 

3.3276C8 

2.947689 

1.406455 

-2.989642 

-1.162037 

-.515246 

.033682 

.615704 

1.444497 

1.684988 

1.765338 

2.149542 

1,907288 

-2.603851 

-1 .028039 

-.348273 

.243226 

.810289 

1.301237 

1.5X6864 

1.8776C5 

2.116383 

1.896195 

-1.975676 

-.760129 

-.158782 

.399574 

.853672 

1.162463 

1.365254 

1.784550 

1.886074 

1.675C9G 

-1.232475 

-.431047 

.051963 

.519137 

.8159C0 

1.057363 

1.218719 

1.431338 

1.416565 

1.303411 

-.449881 

-.047302 

.229487 

•550347 

.764423 

.884001 

.930373 

I.C19C15 

1.C2451C 

1.006941 

.247051 

.271121 

.356211 

.503910 

.645970 

.694601 

.681087 

.675702 

.692556 

.725895 

•749428 

.472464 

.419244 

.393323 

.435654 

.470539 

.493432 

.497904 

.496201 

.501350 

.695369 

.504723 

.399586 

.347474 

.302145 

.295311 

.267939 

.227060 

.244542 

.289372 

.429761 

.441503 

.379793 

.32C990 

.274796 

.245237 

.204232 

.151471 

.138056 

.142303 

.174226 

.107926 

.1905C9 

.202841 

.197818 

.184113 

.182560 

.171867 

.158509 

.145498 





APPENDIX - Continued 


FLEXIBLE SLCPE : 

INCREHEK7S AT 

LOAD POINTS 

DUE TO JIG 

SHAPE 



•0093265 

•00918C2 

•0062798 

.0068154 

•0041496 


.0015830 

-.0004448 

-.0004662 

•C006476 

.0007980 

.0C735C9 

•0064697 

•0053294 


•0039300 

.0010475 

.0002913 

-.0005372 

-.00C7779 

-.0007730 

-.0003592 

•00*2643 


.0043764 

•0033199 

•0021558 

.0005720 

-.0002385 

-.0004512 

-.0030058 

-.0049405 


-.0053916 

•0029295 

•0022681 

.0014482 

.0004625 

-.0007366 

-.0015592 

-.0021212 


-.0056999 

-.0085141 

-.0090986 

.0003423 

.cooiaio 

-.0003467 

-.0012194 

-.0023196 


-.0040738 

-.0057783 

-.0090879 

-.0109831 

-.0113168 

-.0016547 

-.0020683 

-.0027845 


-.0039609 

-.0054903 

-.0077422 

-.0097738 

-.0117936 

-.0159548 

-.0172270 

-.0054048 


-.0059996 

-.0070136 

-.0083952 

-.0104494 

-.0129421 

-.0145660 

-.0162609 

-.0197281 


-.0200938 

-.0110115 

-.0117441 

-.013X881 

-.0149815 

-.0176375 

-.0186922 

-.0202817 


-.0211958 

-.0211937 

-.02106r04 

-.019C899 

-.0198061 

-.0212909 

-.0224324 

-.0234134 


-.0243034 

-.0253232 

-.0282140 

-.0290716 

-.0281724 

-.0237654 

-.0241633 

-.0252019 


-.0255808 

-.0258463 

-.0263332 

-.0269567 

-.C292188 

-.0290195 

-.0293416 

-.0263197 


-.0266395 

-.0269513 

-.0269314 

-.0263350 

-.0264693 

-.0265132 

-.0256829 

-.0258345 


-.0261660 



FLEXIBLE SLCPE 

INCREMENTS AT 

LOAD POINTS 

DUE TC ALPHA 




.2718150 

.2653555 

.23954C7 

.2004067 

.1269749 


•0412762 

-.1042503 

-.1480572 

-.1612471 

-.1637413 

•2184813 

.1925376 

.1595792 


.1195562 

•0493260 

-.0323341 

-.1028083 

-.1515475 

-.1806034 

-.2044741 

.1623259 


.1437793 

.1042767 

.0630070 

-.0123659 

-.C846672 

-.1305597 

-.2383490 

-.3234186 


-.3452587 

.1257243 

•0817493 

.0316665 

-.0216588 

-.0948301 

-.1661607 

-.2091656 


-.3364940 

-.4250158 

-.4429938 

.0511167 

.0077199 

-.0576662 

-.1397261 

-.2073474 


-.2938092 

-.3562284 

-.4431237 

-.4854933 

-.4917673 

-.1076195 

-.1651017 

-.2438254 


-.3176108 

-.3689685 

-.4309977 

-.4856519 

-.5599290 

-.6386993 

-.6645427 

-.3X11987 


-.3683606 

-.4418872 

-.4935266 

-.5475559 

-.6000526 

-.6366479 

-.6923714 

-.74569C8 


-.7546954 

-.5492655 

-.5885708 

-.6365889 

-.6937276 

-.7586663 

-.77X61X8 

-.7955206 


-.8124320 

-.8086727 

-.8043788 

-.7659422 

-.780 0528 

-.8272612 

-.8656580 

-.0997523 


-.9282906 

-.9703982 

-1.1254435 

-1.1704522 

-1.1213711 

-.9742791 

-.9781178 

-1.0145050 

-1.0472356 

-1.0876102 

-1.1276822 

-1.1736419 

-1.2057610 

-1.3458948 

-1.3210345 

-1.1878695 

-1.2101109 

-1.2365246 

-1 .2868476 

-1.3320485 

-1.3877344 

-1.4177012 

-1.3918812 

-1.4161577 

-1.4417986 



FLEXIBLE SLCPE : 

INCREMENTS AT 

LCAD POINTS 

DUE TC DELTA 




-.0002942 

-.0002654 

-.0003076 

-.0002935 

-.CC02632 


-.0012X47 

-.0110868 

-.0173026 

-.0192107 

-.0198353 

-.0002903 

-.0002683 

-.0001929 


-.0000675 

-.0007577 

-.0044752 

-.0106517 

-.0178087 

-.0242425 

-.0298412 

.0002191 


•0010895 

.0011256 

•0002690 

-.0034019 

-.C074713 

-.0114313 

-.0372757 

-.0697415 


-.0814037 

•0005580 

.0007344 

-.0007055 

-.CC34182 

-.0089266 

-.C169C98 

-.0241130 


-.0586673 

-.0955201 

-.1065206 

.0010400 

-.0014968 

-.0056345 

-.0112731 

-.0183115 


-.0356909 

-.0524458 

-.0829878 

-.C990O58 

-.1024402 

-.0047330 

-.0089656 

-.0162095 


-.0277846 

-.0420577 

-.062 9740 

-.0038032 

-.0804850 

-.1336482 

-.1496426 

-.0284944 


-.0341563 

-.0456250 

-.0595292 

-.0773259 

-.0983646 

-.1121435 

-.1078841 

-.1632385 


-.1675837 

-.0690619 

-.0761233 

-.0904911 

-.1027895 

-.1220767 

-.1299337 

-.1416876 


-.1505454 

-.1565484 

-.1566422 

-. 1225493 

-.1247763 

-.1254691 

-.1266107 

-.1248482 


-.1217907 

-.1104183 

-.1056200 

-.1019926 

-.1C63898 

-.1234189 

-.1210463 

-.1169590 


-.1109565 

-.1048641 

-.0980896 

-.0903709 

-.0760489 

-.0711093 

-.0731299 

-.C864065 


-.0052106 

-.0806223 

-.0742764 

-.0685278 

-.C619235 

-.0566827 

-.0572186 

-.053C778 


-.0496569 



FLEXIBLE SLCPE ! 

INCREMENTS AT 

LOAD POINTS 

DUE TO AN, 

G 




-.0268041 

-.0259856 

-.0225443 

-.0100043 

-.01C7137 


-.0035940 

•0056305 

•0C84406 

.0108165 

.0113741 

-.0193221 

-.0168675 

-.0137837 


-.0100608 

-.0042747 

.0012337 

.0055197 

.0062619 

.0103155 

.0126788 

-.0138108 


-.0113385 

-.0083431 

-.0048901 

.0003562 

.0043495 

.0070338 

.0110650 

.0154848 


.0170706 

-.0079015 

-.0056559 

-.0029126 

.0002944 

•0C43538 

.0074153 

.0093750 


.0137153 

.0175523 

•0186144 

-.0016061 

.0002251 

.0027609 

.0050058 

.0081602 


•0112828 

.0132437 

.0159840 

.0172588 

•0171361 

•0046172 

.0063574 

.CC86204 


•0108256 

.012493Q 

•0142663 

.0153492 

.0161658 

.0161123 

.0159463 

.0104969 


.0115695 

.0130912 

•0142618 

.0151451 

.0157920 

•0160051 

.0157588 

.0154953 


•0154811 

.0149438 

.0152722 

.0154900 

.0157221 

.0158398 

.0155771 

.0153745 


.0151743 

•0151365 

•0151684 

.0154060 

.0154449 

.0152302 

•0149953 

.0146808 


.0143669 

•0139186 

.0115150 

.01C818B 

.0115696 

.0140530 

.0138967 

.0135250 


.0130663 

•0125265 

.0120453 

.0114763 

.0094645 

•0088746 

•C0925I9 

.0110281 


.0109195 

•0106546 

•0100214 

.0092955 

•0084688 

.0080720 

.0085677 

.0082501 


.0078779 




APPENDIX - Continued 


LCAD FACTOR = .99375 



TRIM VALLES i 

2 

ITEBATICNS) 



CA = 

.004759 

ALPHA 

= 

•046496 

PHI * 

c.occooc 

cn - 

.C79145 

DELTA 


-.001151 

ANf 6 = 

.992585 

CM = 

C. 000000 

THETA 

= 

.048496 

GC/2V = 

.000000 


DERIVATIVE CCMRIBUTICNS 



JIG SHAPE 

ALPHA 

DELTA 

QC/2V 

AN. G 

QDOT 

PACH 

OYN PRESS 


CA 


-.0396468 

.0012515 

-.CC99517 

.0001567 

-.0012869 

.0008278 

-.C00CCC47 

CN 

.014051 

1.119244 

.C2C705 

.194543 

.010917 

-.032467 

-.011269 

-.0000255 


CM 

.005452 

.C43187 

-.016252 

-.372356 

-.0C7624 

.001052 

-.0CC896 

.000C17^ 



STABILITY 

CEP IVATIVES 









U 

LCCT 

ALPHA 

ALPHA COT 

THETA 

0 

QDOT 

DELTA 

ALTITUDE 

CA 

.0016779 

-.0006172 

-.0396468 

-.569C303 

-.CCCC076 

.5790866 

-.0012869 

.0012515 

1.32361E-O0 

CN 

-.060885 

-.042995 

1.119244 

-41.03C6C6 

-.000526 

41 .225150 

-.032467 

.020705 

7.23373E-07 

CM 

.013285 

.C30C27 

.043187 

28.655362 

.000367 

-29.027718 

.001052 

-.018252 - 

4.91711E-07 

X 

-.0C1728 

.CC0Q95 

.0C612O 

.0C1963 

.COOCOl 

-.001930 

.0C0199 

-.000193 

3.28540E-08 

z 

-.015037 

.0C6637 

-.172782 

.136749 

.OOOC01 

-.137397 

.0C5012 ■ 

-.003196 

4.68681E-07 

M 

.287332 

.471844 

.678637 

9.721432 

.0C5772 

-9.847756 

.016531 

-.286803 - 

7.72668E-06 


STATIC FAPARETEPS 


CM ALPHA/CN ALPHA 

- 

.026566 

STATIC MARGIN 

- 

-.091774 

MANEUVER MARGIN 

= 

-.063460 

DELTA/U 

= 

.762541 

DBLTA/AN, G 

= 

-.265611 


TOTAL 

TRIM CCNCITION SLOPES 

AT SLOPE POUTS 




-.0020586 

.0217548 

.0650027 

.0721600 

.0677658 

.0584700 

.0391715 

.0218C80 

-.0198535 

-.0640676 

-.0173102 

.0054710 

.0166357 

.0223172 

.0250957 

.0253133 

.0236225 

.0206882 

.0164255 

.0115561 

-.C4C7771 

-.0091254 

-.0025005 

.0047652 

.0093337 

.0117247 

•0126286 

.0129417 

.0118786 

.0090395 

-.0524002 

-.0241931 

-.0143279 

-.CC61938 

-.0008480 

.0026211 

.CC542E5 

.0065721 

.0070041 

.0073536 

-.0647082 

-.0367440 

-.0204247 

-.0143891 

-.0089509 

-.0049122 

-.0014923 

.0009690 

.0027601 

.CC42169 

-.0766635 

-.0474647 

-.0314577 

-.0236140 

-.0175042 

-.0125537 

-.0090044 

-.0056184 

-.0034212 

-.0016896 

-.0850557 

-.0638657 

-.0417912 

-.0301311 

-.0240834 

-.C207QOC 

-.0174310 

-.0138160 

-.0119256 

-.0101637 

-.0978100 

-.0723314 

-.0569273 

-.0477079 

-.0418801 

-.0354738 

-.0309467 

-.0292457 

-.0272503 

-.0233828 

-.0665331 

-.0612836 

-.C5661C0 

-.0504553 

-.0461084 

-.0429040 

-.0394450 

-.0369537 

-.0351936 

-.0336691 

-.0209370 

-.0400142 

-.0420570 

-.0427407 

-.0432207 

-.0435807 

-.0439389 

-.0441645 

-.0442936 

-.0442990 

.0066955 

-.0252318 

-.0367012 

-.0433458 

-.0433476 

-.0433460 

-.0433572 

-.0433594 

-.0433567 

-.0433567 




TOTAL 

TRIM CCKC 

ITICN SLOPES 

AT LCAC POINTS 




-.0436302 

.1331363 

.0554884 

.0755297 

.0738253 

.0657601 

.0510246 

.0254600 

.0046302 

-.0542782 

-.0404576 

-.C031C56 

.0133686 

.0200093 

.0247193 

.0260647 

.0248239 

.C225010 

•0186945 

.0137691 

-.0615112 

-.0174440 

-.0048862 

.0023790 

.0079593 

.0110990 

.0127490 

.0132884 

.C1265C6 

•0105662 

-.0688587 

-.0334707 

-.0183161 

-.C09C0C0 

-.0021438 

.0013085 

.0050326 

.0061520 

.0068567 

.0071206 

-.0866681 

-.0458635 

-.0260509 

-.0170926 

-.0105223 

-.0064905 

-.0025338 

.0001094 

.0022373 

.0036779 

-.0961360 

-.0585503 

-.0374574 

-.0259980 

-.0200820 

-.0139061 

-.0104473 

-.CC68051 

-.0042518 

-.0024812 

-.0882575 

-.0722216 

-.0504121 

-.0334187 

-.0254524 

-.0217348 

-.01 87572 

-.0151476 

-.0125261 

-.0107558 

-.0042833 

-.0371175 

-.0611572 

-.0485862 

-.0432940 

-.0368230 

-.0305567 

-.0283125 

-.0274670 

-.0235912 

-.0635242 

-.C648090 

-.0599264 

-.0539971 

-.0477563 

-.0450994 

-.0411916 

-.0380085 

-.0359937 

-.0340499 

-.0173197 

-.0366564 

-.0417364 

-.0425208 

-.0431453 

-.0434177 

-.0438509 

- .044C692 

-.0442386 

-.0442742 

.0256397 

-.0114377 

-.0340408 

-.0434761 

-.0433465 

-.0433451 

-.0433468 

-.0433562 

-.0433570 

-.0433554 
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Program Listing 

The source program listing is presented on the following pages. The listing for 
subroutine MATRIX is not given, but its fxmction can be performed by any standard rou- 
tines for matrix multiplication and inversion. In the listing presented, any CALL MATRIX 
statement having a first argument of 20 calls for multiplication of the matrix starting at 
the address named in the fifth argument into the matrix starting at the address named in 
the seventh argument, the results to be stored starting at the address named in the ninth 
argument. Any CALL MATRIX statement having a first argument of 10 calls for inversion 
of the matrix starting at the address named in the fifth argument, the inverted matrix to 
be stored starting at the same address, and the determinant to be stored at the address 
named in the seventh argument. The determinant is not used in the present program. 

The remaining arguments transmit instructions and dimensions required by subroutine 
MATRIX. 


65 



onooonooonnonoooooonoooonnonononoooooononnnnnoooo 


APPENDIX - Continued 


PROGRAM LISTING 


HAS DIMENSIONS M*M AND IS OF THE FORM SO THAT S(ROW,COU 
♦FORCE(COL) = SLOPEIROWI* TAPE MUST BE WRITTEN WITH COLUMN 
INDEX ADVANCING MOST RAPIDLY AND NO END FILE CODE BETWEEN 
MATRICES. 


PROGRAM AEOERIVC INPUT, OUTPUT tTAPE5=INPUT, TAPE6=aUTPUT»TAPE 1 tTAPEZ I AEO 

AED 

TAPEl MUST CONTAIN THREE AERO IC MATRICES FOR MACH = AM, MACH = AEO 
AM+OM, ANO MACH * AM-CM. EACH MATRIX HAS DIMENSIONS M^M AND IS AED 
OF THE FORM SO THAT A(ROW,COU*SLOPE ( COL ) = FORCE(ROWl/Q. TAPE AED 
MUST BE WRITTEN WITH COLUMN INDEX ADVANCING MOST RAPIDLY ANO NO AEO 
END FILE COOES BETWEEN MATRICES. AEO 

AED 

TAPE2 MUST CONTAIN TWO STRUCTURAL SLOPE IC MATRICES FOR SLOPES AT AEO 
PANEL SLOPE POINTS, ANO SLOPES AT PANEL LOAD POINTS. EACH MATRIXAEO 

AEO 
AED 
AED 
AED 
AEO 
AEO 
AED 
AED 
AEO 
AEO 
AEO 
AEO 
AEO 
AEO 
AEO 
AEO 
AED 
AED 
AED 
AED 
AEO 
AED 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 


IFLEX = 0 FOR RIGID CASE (REQUEST TAPEl ONLY I 
IFLEX = 1 FOR ELASTIC CASE (REQUEST TAPEl AND TAPE2) 

LIST = 0 DELETES LISTING BELOW 

LIST = 1 LISTS ELASTIC AIRLOADS ANO FLEXIBLE SLOPE INCREMENTS 
DUE TO JIG SHAPE, ALPHA, DELTA AND AN, G PLUS TOTAL TRIM 
CCNDITION SLOPES 

LFT = 0 TRIMS WITH WINGS LEVEL AT THE INSTANT THAT FLIGHT PATH 
IS LEVEL 

LFT = 1 TRIMS IN CONSTANT ALTITUDE BANKED TURN IF LOAD FACTOR IS 
GREATER THAN G/32.174 
NC = NUMBER OF CG LOCATIONS 

NW = NUMBER OF AIRPLANE WEIGHTS FOR EACH CG 
NQ = NUMBER OF Q VALUES FOR EACH WEIGHT 
NN = NUMBER OF LOAD FACTORS FOR EACH 0 

CAF = INCREMENT IN AXIAL FORCE COEFFICIENT (SKIN FRICTION, ETC) 

TO BE ADDED TO THAT CALCULATED DUE TO SURFACE PRESSURES 
XCG = X COORDINATE OF CG (FROM NOMINAL ORIGIN, X POSITIVE FORWARO)AED 330 

AED 340 

INPUT ARRAYS AED 350 

TSR(H,1) = SLOPES AT SLOPE POINTS DUE TO TWIST ANO CAMBER AEO 360 

TSR(M,3) = SLOPES AT SLOPE POINTS DUE TO UNIT DELTA (RAO) AED 370 

TLRJ(M) = SLOPES AT LOAD POINTS DUE TO TWIST AND CAMBER AEO 380 

TLRC(M) = SLOPES AT LCAO POINTS DUE TO UNIT DELTA (RAD) AED 390 

XSC(M) = X COORDINATE OF SLOPE POINTS (FROM NOMINAL ORIGIN, AEO 400 

X POSITIVE FORWARD) AEO 410 

XLC(M) = X COORDINATE OF LOAD POINTS AED 420 

FI(M,1) = WEIGHTS ACTING AT PANEL LOAD POINTS AEO 430 

FL(NN) = LOAD FACTOR AEO 440 

LOAD FACTOR IS NORMAL TO FLIGHT PATH, IN SEA LEVEL G UNITS AEO 450 

AN, G IS NORMAL TO THE XY PLANE OF BODY AXES, IN SEA LEVEL G UNITSAEO 460 

BODY AXES ARE FIXED TO A MATERIAL POINT AT THE REFERENCE POINT AEO 470 

OF THE STRUCTURAL SLOPE MATRIX AED 480 

AED 490 
AED 500 

COMMCN A(120,120) ,S( 120 , 12C ) , TSR( 1 2C , 7 ) , FR A( 120 , 7 ) , XS3( 120 ) , XLO( 1 2AED 510 
10 ) ,XS(120) ,XL( 120),FI ( 120 , 2 ) , FL ( 2 5 ) , F E A ( 12C , 9 ) , F I A ( 120 , 3 ) , CN( 8 ) , CMAE 0 520 


(8), 

TLE 

(120, 

9) ,CA(8) 

,TLRJ(120),TLR0( I20),SCA(9) 

,SCN(9),SCM(9) 

,SX(9AE0 

530 

\)fSl 

(S) 

,SM(9 

) 



AED 

540 

REAL 

(5 

,54 ) 

IFLEX, LI 

ST,LFT,M,NC,AM,DM, SREF,CBAR 

,CAF 

AEO 

550 

DO 1 

K= 

1,3,2 




AEO 

560 

READ 

(5 

,55) 

(TSR( I,K 

), I=1,M) 


AEO 

570 

READ 

(5 

,55) 

(TLRJ( I) 

,I=1,M) 


AEO 

580 
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APPENDIX - Continued 


RE^C {5f55) (TLRO( n f 
REhlbiO 1 

READ (1) ( (A(I»J> 

IR=1 

DO 2 1 = 1, M 
TSR( I,2)*l. 

CGHPUTE RIGID LOADS DUE TO JIG SHAPE, ALPHA AND DELTA 
DO 3 K»l,3 

CALL MATRIX ( 20 ,M ,M, 1 , A, 120 , TSR ( 1 , Kl , 120,FRA( 1 , K ) , 120 I 
CONTINUE 

READ (5,551 ( XSO ( 1) , 1 = 1 , M ) 

READ (5,55) ( XLC ( I),I = 1,M) 

SRC2=SREF/2* 

DO 53 IC=l,NC 

COMPUTE RIGID LOADS DUE TO QC/2V 

READ (5,56) NW,XCG 

DO 4 1=1, M 

XS( I ) = XS0( I)-XCG 

XL( I ) = XL0( n-XCG 

TSR( I,4)=-XS( I I *2. /C BAR 

IF ( IR.EQ.l) GO TO 5 

REWIND 1 

READ (1) ( ( A( I , J) , J=1,M) ,I = 1,M) 

IR = 1 

CALL MATRIX ( 2C , M , M, 1 , A , 1 20 , T S R ( 1 , 4 ) , 1 20 , FRA ( 1 , 4 ) , 1 20 ) 

DO 52 IW=1,NW 

READ (5,56) NQ , W , Y I , Y I , Y2 

ESTABLISH INERTIAL LOADS DUE TO LOAD FACTOR AND ODOT 

IF ( IFLEX.EQ.O ) GO TO 7 
READ (5,57) (FI( I!,l) ,I=1,M) 

DO 6 1=1, M 
FI(I,1)=-FI(I,1) 

FI ( I ,2)=FI (I ,l )*XL( I )/32.174 
DO 51 IQ=1,NQ 

READ (5,56) NN ,0 , V ,G , RHC , RHOH , AH 
READ (5,55) (FL( I ),I = 1,NN) 

WRITE (6,58) AM,XCG,W,Q 

READ BASIC A AND S MATRICES 

IF ( IR.EQ.l) GO TO 8 
REWIND 1 

READ (1) .((A(I,J),J=1,M),I = 1,M) 

IF ( IFLEX.EQ.O) GO TO 13 
REWIND 2 

READ (2) ((S(I,J) ,J=I,M),I = 1,M) 

FIND AIR LOADS DUE TO INERTIAL DISTORTION 

DO 9 KI=1,2 
K»Kl+4 

CALL MATRIX ( 20 , H , M, 1 , S , 120, F I ( 1 , K I ) , 1 2C , TSR ( 1 , K ) , 12C ) 
CALL MATRIX ( 20, M ,M, 1 , A ,120 ,TSR( 1 , K) , 120 ,FRA ( 1 , K) , 120 ) 

COMPUTE BASIC AERGELASTIC CORRECTION MATRIX 


AED 590 
AEO 600 
AED 610 
AEO 620 
AEO 630 
AEO 640 
AEO 650 
AED 660 
AED 670 
AEO 680 
AED 690 
AED 700 
AED 710 
AED 720 
AED 730 
AEO 740 
AEO 750 
AED 760 
AED 770 
AEO 780 
AED 790 
AEO 800 
AED 810 
AEO 820 
AEO 830 
AEO 840 
AEO 850 
AEO 860 
AEO 870 
AEO 880 
AEO 890 
AEO 900 
AEO 910 
AEO 920 
AEO 930 
AED 940 
AED 950 
AED 960 
AED 970 
AED 980 
AEO 990 
AEDIOOO 
AEOlOlO 
AE01020 
AED1030 
AE01040 
AE01050 
AED1060 
AE01070 
AED1080 
AE01090 
AEOllOO 
AEDlllO 
AED1120 
AE01130 
AE01140 
AE01150 
AED1160 
AE01170 
AE01180 
AE01190 
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c 



AED1200 


CALL MATRIX ( 20 , M ,M , M ♦ A , 120 t S , 120 » A, 1 20 ) 


AED1210 


DO 10 J=lfM 


AE01220 


00 1C 1=1, M 


AE01230 


A(I,J)=-Q*A(I,J) 


AED1240 


IF (KEC.J) A( I, J)=A( I,J)4^K 


AE01250 

10 

CONTINUE 


AE01260 


CALL MATRIX ( 10 , M ,M , 2 , A , 120 , DETERM ) 


AE01270 

c 



AE01280 

c 

APPLY AEROELASTIC CORRECTION TO FIRST SIX SETS OF AIR 

LOADS 

AED1290 

c 



AE01300 


DO 11 K=l,6 


AE01310 

11 

CALL MATRIX ( 20 , M ,M , 1 , A , 1 20 , FRA ( 1 , K > , 1 2C , F E A ( 1 , K ) , 1 20 ) 


AED1320 

c 



AED1330 

c 

SUM AIR AND INERTIAL LOADS FOR INERTIAL LOAD SOURCES 


AE01340 

c 



AED1350 


DO 12 KI=1,2 


AED1360 


K=Kl+4 


AED1370 


DO 12 1=1, M 


AE01380 

12 

FIA(I,KII=FEA( I,K)+FI( I,KI)/Q 


AE01390 


Kl = 6 


AE01400 


GO TO 16 


AE01410 

13 

DO 15 1=1, M 


AED1420 


DO 14 K=5,6 


AED1430 


TSR( I,K)=0, 


AE01440 

14 

FRA( I,K)=0. 


AE01450 


00 15 K=l,6 


AE01460 

15 

FEA(I,K)=FRA(I,K) 


AED1470 


CN(5)=C. 


AE01480 


CN(6)=0. 


AE01490 


CM(5)=0. 


AE01500 


CM(6)=0. 


AE01510 


Kl = 4 


AE0152D 

16 

II 

00 

o 

o 


AE01530 

C 



AE01540 

C 

COMPUTE CN AND CM DUE TO FIRST SIX LOAD SOURCES 


AED1550 

C 



AE01560 


CN(K)=C. 


AE01570 


CM(K)=0. 


AE01580 


DO 17 1=1, M 


AED1590 


CN(K)=CN(K)+FEA(E,KI 


AE01600 

17 

CM(K)=CM(KI + FEA( I,K)*XL( I» 


AE01610 


CN(K)=CNCK>/SR02 


AED1620 

18 

CM(K)=CM(K)/(SR02*CBARI 


AED1630 


IR=0 


AED1640 


IF (LIST.EQ.O) GO TO 19 


AE01650 


WRITE (6,59) (FEA(I,1),I=1,M) 


AED1660 


WRITE (6,60) (FEAI 1,2) ,I = 1,M) 


AED1670 


WRITE (6,61) (FEA(I,3) ,I = 1,M) 


AE01680 


IF (IFLEX.NE.O) WRITE (6,62) ( FEA ( I , 5 ) , I =1 , M ) 


AED1690 

19 

DO 50 IN=1,NN 


AED1700 

C 



AED1710 

C 

COMPUTE TRIM CONDITIONS 


AE01720 

C 



AE01730 


CALL TRIM (CN,CM,H,FL( IN),LFT,Q,SREF,CBAR,V,G,YI,Y1,Y2 

,AT,DT,ANT 

,CAE01740 


1NT,CMT,QT,PR,PHI ,T,I TER) 


AED1750 


IF (ITER. GE. 100) GO TO 49 


AED1760 

c 



AED1770 

c 

COMPUTE TRIM SETS OF RIGID SLOPES. RIGID AIR LOADS AND 

ELASTIC 

AED1780 

c 

AIR LOADS 


AE01790 

c 



AED1800 
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20 


21 

C 

C 

C 


22 


C 

C 

C 


23 


C 

C 

C 

24 

25 

26 
C 
C 
C 
C 

27 


28 


00 20 I=1,M 

TSR( I,7)=TSRC Ifl)+AT+0T*TSR( If 3)+QT*TSR( I,4J+ANT^TSRC If 5)-PR*TSR( 

1 f 6 ) 

FRA( i,7)=FRA< I»1 >+AT*FRAC I,2)+DT*FRA( I,3)+QT*FRA( 1 , 4 > +ANT*FRA( I» 5 
1-PR*FRAII,6) 

F6A(If9)=:FEA(I,l)+AT*FeA(I,2l+0T*FEA( I,3) + QT*FEA( 1 , 4) +ANT*FEA( 1,5 
1-PR*FEA(I,6I 
IF (IFLEX.EQ.O) GC TO 25 
DO 21 1=1, M 

FIA(I,3)=FEA(I,9)+(ANT*FI( I,1)-PR*FIC I , 2 ) ) /Q 

CORRECT TRIM RIGID AIRLOADS FOR AEROELASTIC I TY AT HIGHER Q 

QP=1.05*Q 
REhINC 1 

READ U) ((AUfJ), JslfP)fI=l,M) 

CALL MATRIX < 20 , M, M, M, A, 120, S, 120 , A, 1 20 ) 

DO 22 J=l,M 

DO 22 1=1, M 

A( I, J)=-QP*A(It 

IF (I»EQ.J) A(I,JI=A(I,J)>1. 

CONTINUE 

CALL MATRIX 110 , M , M, 2 , A , 120 , OETERM ) 

CALL MATRIX ( 20 , M , M , I , A , 120 , FRA ( 1 , 7 1 , 1 20 , F EA ( 1 , 8 I , 1 20 ) 

CORRECT TRIM RIGtC AIRLOADS FOR AEROELASTIC I TY AT LOWER Q 

OP=,95*Q 
REWIND I 

READ m nA(l,J>,J=l,M),I = l,M) 

CALL MATRIX ( 20 , M ,M , M , A , I 20 , S , 120 , A, 1 20 ) 

00 23 J=1,M 
DO 23 1=1, M 
A(I,J)*-QP*A(I,J) 

IF (I.EC.J) A(I, J)=A( IfJI^l. 

CONTINUE 

CALL MATRIX Cl 0 , M ,M, 2 , A , 1 20 , OETERM ) 

CALL MATRIX ( 20 , M , M , I , A , 1 20 , FR A ( 1 , 7 ) , 1 20 , F E A ( 1 , 7 ) , 1 20 ) 

COMPUTE ELASTIC AIRLOADS DUE TO Q 
DO 24 1=1, M 

FEAC I,8> = CFEA( I,8>-FEA( 1,7) )/(.l^Q) 

GO TO 27 
DO 26 1=1, M 
FEA( I,8)=0. 

USING TRIM RIGID SLOPES, COMPUTE RIGID AND ELASTIC AIR LOADS 
AT HIGHER MACH 

READ (1) ( ( A(I , J) , J=1,M) , I=1,M) 

CALL MATRIX (20,M,M, 1 , A , 120 , TSR ( 1 , 7 ) , 1 20 , F £ A ( 1 , 7 ) , 1 2<^ ) 

IF ( IFLEX.EQ.O) GO TO 29 

CALL MATRIX C 20 , M , M, M , A , 1 20 , S , 1 2C , A, 1 20 ) 

DO 28 J=l,M 
DO 28 1 = 1, H 
A(I,J)=-Q*A(l,J) 

IF (l.EQ.J) A(l,J)=A(I,J)+l« 

CONTINUE 

CALL MATRIX ( 10 , M ,M , 2 , A , 120 , DETERM ) 

CALL MATRIX C 20 , M , M , 1 , A , 120 , FE A ( 1 , 7) , 1 20 ,FE A ( 1 , 7 ) , 120 ) 


AE01810 

IAED1820 

AE01830 

)AE01840 

AE01850 

)AEDI860 

AED1870 

AE01880 

AE01890 

AE01900 

AE01910 

AE01920 

AE01930 

AE01940 

AED1950 

AED1960 

AE01970 

AE01980 

AE01990 

AE02000 

AE02010 

AED2020 

AED2030 

AED2040 

AE0205D 

AE02060 

AED2070 

AE02080 

AED2090 

AE02100 

AED2110 

AED2120 

AE02130 

AE02140 

AE02150 

AED2160 

AED2170 

AE02180 

AE0219O 

AE02200 

AED2210 

AED2220 

AED2230 

AED2240 

AED2250 

AED2260 

AED2270 

AE0228^ 

AE02290 

AED2300 

AE02310 

AE02320 

AED2330 

AED2340 

AED2350 

AED2360 

AE02370 

AE02380 

AE02390 

AED2400 

AE02410 
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c 


AE02420 

c 

USING TRIM RIGID SLOPES. COMPUTE RIGID AND ELASTIC AIR LOADS 

AE02430 

c 

AT LOWER MACH 

AE02440 

c 


AED2450 

29 

READ (11 ((A(I«JI.J:=1«M).1»1«M| 

AE02460 


CALL MATRIX (20.MtH. 1. A«120tTSR(1.7) .120. TLE( 1. 1 ) .120 ) 

AE02470 


IF (IFL£X*EQ.0I GO TO 31 

AE02460 


CALL MATRIX (aO.H.M.M. A.120.S.120. A. 1201 

AED2490 


00 30 

AE02500 


00 30 1«1«M 

AE02510 


A( I. J)=-Q*A(I. Jl 

AE02520 


IF n.EG.J) AU« JMA( 1.J|4>1. 

AE02530 

30 

CONTINUE 

AE02540 


CALL MATRIX (lO.M «M.2. A . 120. OETERH ) 

AED2550 


CALL MATRIX (20.M.M. 1. A. 120.TLE( 1 . II . 120.TLE( I. 1 1 • 120 1 

AE02560 

C 


AE02570 

C 

COMPUTE ELASTIC AIRLOAOS DUE TO MACH 

AEO250O 

C 


AEO2590 

31 

00 32 I=l»M 

AED2600 

32 

FEAC I,7)=(FEACI.7I-TLE( I,in/(2.*0M) 

AE02610 

C 


AE02620 

C 

FIND CN AND CM DUE TO C AND MACH 

AED2630 

C 


AED2640 


00 34 K=7,8 

AED2650 


CN(K)=0. 

AE02660 


CM(K)=C. 

AE02670 


00 33 1 = 1, M 

AE02680 


CN(K|=CN(K)+FEA( I.K) 

AED2690 

33 

CM(K»=CM(K)+FEA( I,K|^XL( I ) 

AED2700 


CN(KI=CN(K)/SR02 

AED2710 

34 

CM(KI=CM(K)/(SR02*CBAR) 

AE02720 


IF (IFLEX.EQ.O) GO TO 36 

AED2730 

C 


AE02740 

C 

COMPUTE TRIM CONDITION SLOPES AT SLOPE POINTS 

AED2750 

C 


AED2760 


CALL MATRIX (20 .M ,M, 1 , S . 120 , F lA ( 1 . 3 ) . 120 ,TLE ( 1 , 9 1 , 120 I 

AE02770 


DO 35 1=1, M 

AED2780 

35 

TSR(I,7)=TLE( I ,9I*Q+TSR(I,1)+0T*T5R( 1,31 

AE02790 


GO TO 39 

AE02800 

36 

DO 37 1 = 1, M 

AE02810 

37 

TSR( I,7) = TSR( I ,l)+DT*TSR(I ,3) 

AE02820 


00 38 K=l,9 

AED203O 


00 38 1=1, M 

AE02840 

38 

TLE( I.K MO. 

AED2850 


GO TO 43 

AE02860 

C 


AE02870 

C 

COMPUTE SLOPES AT LOAD POINTS FOR ALL LOAD SOURCES 

AE02880 

c 


AE02890 

39 

READ (21 ((S(I.J) tJ=l,MI,Ml,M) 

AE02900 


DO 40 K=l,8 

AE02910 


IF (K.EQ.5.0R.K.EQ.6) GO TO 40 

AED2920 


CALL MATRIX ( 20 .M ,M, 1 , S . 120 , FEA ( 1 , K) , 1 20. TLE ( 1 , K ) , 120 1 

AED2930 

40 

CONTINUE 

AED2940 


DO 41 K=5,7 

AED2950 


Kl=K-4 

AED2960 


L = K 

AE02970 


IF (K,EQ.7I L=9 

AE02980 

41 

CALL MATRIX ( 20 ,M ,M , 1 , S , 120 , F I A ( 1 , K1 ) , 1 20 , TLE ( 1 , L ) , 120 1 

AE02990 


CALL MATRIX ( 20 ,M , M , 1 , S , 120 ,FE A (1 , 9) , 1 20 ,F I A ( 1 , 3 1 , 1 20 I 

AE03000 


DO 42 K=l,9 

AE0301C 


00 42 1=1, M 

AED3020 
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TLEnfK)^TLEntKl*g 
IF <K.NE.8) GO TO 42 
TLEC I.8)*TLen,8)+FIA( I»3) 

42 CONTINUE 

IF (IFLEX.EQ.0.OR.LIST.EQ*0) GO TO 43 

IF IIN.NE.n GO TO 43 

WRITE (6t63) (TLE ( 1 1 1 ) 1 1^=1 1 M) 

WRITE (6,64) (TLE(I,2I f I=I,H) 

WRITE (6,65) (TLE(I,3) ,I»1,M) 

WRITE (6,66) (TLE( 1, 5) ,I^1,M) 

C 

C COMPUTE AXIAL FORCE COEFFICIENT AND DERIVATIVES 

C 

43 CAT=0. 

DO 44 1=1, M 

TLE( It3) = TLE(I,3)+TLRO(I) 

TLE( I,9) = TLE( I,9)>TLRJ( I ) + DT*TLRD( I) 

44 CAT=CAT+(TLE( I,9)*FEA( 1,9) )/SR02 
00 46 K=2f8 

CA(K)=0. 

DO 45 1=1, M 

45 CA(K)=CA(K)+(TLE( I ,9 ) *FE A ( I , K ) +TLE ( I , K ) *FE A ( 1,9) )/SR02 

46 CONTINUE 
CAT=CAT+CAF 

WRITE (6,67) FL( IN) 

WRITE (6,68) ITER,CAT, AT,PHI,CNT,DT, ANT,CMT,T,QT 
WRITE (6,69) (CA(K),K=2f8) ,CN,CM 

C 

C COMPUTE COMBINED STABILITY DERIVATIVES 

C 

FN=2.*V*V/(32.174*CBAR) 

VOG=V/32.174 

CSA=CCS(AT) 

SNA=SIN(AT) 

ST=SIN(T) 

CP=C0S(PHI ) 

SCA(1)=FN*QT*CSA+CA(5)+AM*CA(7)+2.*Q*CA(8) 

SCN( 1)=FN*QT*CSA*CN( 5)+AM*CN(7)+2.*Q*CN(8) 
SCM(1)=FN*QT*CSA*CM(5)4AM*CM(7)+2.*Q*CM(8) 
SCA(2)=-VOG*SNA*CA(5) 

SCN( 2)=-V0G*SNA*CN(5) 

SCM(2)=-V0G*SNA*CM(5 ) 

SCA( 3)=CA( 2)-FN*QT*SNA*CA( 5) 
SCN(3)=CN(2)-FN*QT*SNA+CN(5) 
SCM(3)=CM(2)-FN*QT*SNA*CM(5) 

SCA(4)=-FN*CSA*CA (5) 

SCN(4)=-FN*CSA*CN(5) 

SCM (4)=-FN*CSA*CM(5l 

$CA(5)=-ST*CP*CA(5)*G/32.174 

SCN(5)=-ST*CP*CN(5)*G/32.174 

SCM(5)=-ST^CP*CM(5)*G/32.174 

SCA(6)=CA(4)+FN*CSA*CA(5) 

SCN(6)=CN(4)+FN*CSA»CN(5) 

SCM(6)=CM(4)+FN*CSA*CM(5) 

SCA(7)=CA(6) 

SCN(7)=CN(6) 

SCM(7)*CM(6) 

SCA(8)=CA(3) 

SCN(B)*CN(3) 

SCM(8)=CH(3) 

SCA(9)=Q*RH0H*CA(8)-AM*AH*CA(7) 


AE03030 

AE03040 

AE03050 

AE03060 

AE03070 

AE03080 

AED3090 

AE03100 

AED3110 

AED3120 

AED3130 

AE03140 

AED3150 

AE03I60 

AED3170 

AED3180 

AED3190 

AE03200 

AE03210 

AE03Z20 

AED3230 

AE03240 

AE03250 

AE03260 

AE03270 

AE03280 

AE03290 

AE03300 

AED3310 

AED3320 

AE03330 

AED3340 

AED3350 

AE03360 

AE03370 

AE03380 

AE03390 

AE03400 

AED3410 

AE03420 

AE03430 

AE03440 

AED3450 

AED3460 

AE03470 

AED3480 

AE03490 

AE03500 

AED3510 

AE03520 

AEO3530 

AED3540 

AE03550 

AE03560 

AE03570 

AE03580 

AE03590 

AED3600 

AE03610 

AED3620 

AE03630 
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APPENDIX - Continued 



SCN 1 9 1 = Q*RHOH*CN ( 8 AM* AH*CN ( 7 ) 

AE03640 


SCF(9l=g*RH0H*CM(8l-AM*AH*CM( 7) 

AE03650 


QM=32,174*Q*SREF/ (W*V) 

AE03660 


QI=Q*SREF*CBAR/YI 

AED3670 


SX ( 1 (SCA ( 1 ) +2. *C AT ) 

AE03680 


SZ C n =-QM* ( SCN( 1 ) *2. *CNT) 

AE03690 


SMa)=QI*(SCM<ll + 2.*CMT) 

AE03700 


DO 47 K=2f8 

AED371C 


SX(K>=-CM*SCACK) 

AE03720 


SZ(K)--CM*SCN(KI 

AE03730 

A7 

SMIK)=QI*SCM|K» 

AE03740 


00 48 K=4,6,2 

AED3750 


SX(K)=SX(K)*C8AR/(2.*V) 

AE03760 


SZ(K)=SZ(K)*CBAR/ (2.*V) 

AE03770 

48 

SM(K)=SM(KJ*CBAR/f2.*V) 

AE03780 


SX(9)=-CM*(SCA(9)+RH0H*CAT) 

AE03790 


SZ ( 9 ) =-CM* ( SCN ( 9 ) +RHOH+CNT > 

AED3800 


SM(9)=gi*($CM(9)*RH0H*CMTI 

AED38in 


WRITE (6,70) SCA,SCN,SCM,SX,SZ,SM 

AED3820 

C 


AE03830 

C 

COMPUTE STATIC PARAMETERS 

AE03840 

c 


AED3850 


CMCN=SCM(3)/SCN(3) 

AE03860 


CMNU=CMCN*ll.<-SCNm / (2.*CNT) )-SCM(l)/(2.*CNT) 

AED3870 


B=32.174*RH0*SREF*CBAR/(4.*W) 

AE03880 


CMNN=CMCN*(1.-B*SCN(6) ) *6*501(6) 

AE03890 


DU=2 •*CNT*CMNU/ ( SCM( 8 ) -CMCN*SCN( 8 ) ) 

AED3900 


DA^-W*CHNN/(Q*SREF*(SCM(8)-CMCN*SCN(8) ) ) 

AE03910 


WRITE (6,71) CMCN,CMNU,CMNN,OU,DA 

AE03920 


IF (LIST. EG. C) GO TO 50 

AE03930 


WRITE (6,72) (TSR(I,7) ,I = 1,M) 

AE03940 


WRITE (6,73) (TLE(I,9) ,I^1,M) 

AE03950 


GO TO 50 

AED3960 

49 

WRITE (6,74) CN,CM 

AE03970 

50 

CONTINUE 

AED3980 

51 

CONTINUE 

AED3990 

52 

CONTINUE 

AED4000 

53 

CONTINUE 

AE04010 


STOP 

AE04020 

C 


AE04030 

54 

FORMAT (3I1,2I3,5F10.0) 

AED4040 

55 

FORMAT (8F10*0) 

AE04050 

56 

FORMAT n3,6F10«0) 

AED4060 

57 

FORMAT (6F12.0) 

AE04070 

58 

FORMAT (1HI,5X,27FAER0ELASTIC DERIVATIVES, M=,F6.3,6H, 

XCG=,F8.3,4AE04080 


IHt W»>FS.0.7H, QB«R-,F8.2I AEDA090 

59 FORMAT I//10X, 39HELAST IC AIRLOADS PER 0 DUE TO JIG SHAPE/t IX,8F12.AEDA100 


1611 AE04110 

60 FORMAT I //lOX, 35HELAST IC AIRLOADS PER Q DUE TO ALPHA/ (1 X, 8F12. 6H AE0A120 

61 FORMAT I//10X,35HELASTIC AIRLOADS PER Q DUE TO DELT A/ 1 IX, 8F12. 6)1 AE0A130 

62 FORMAT I//10X, 35HELASTIC AIRLOADS PER 0 DUE TO AN, G/I IX,8F12.6) ) AEOAIAO 

63 FORMAT I//10X,57HFLEX I 8LE SLOPE INCREMENTS AT LOAD POINTS DUE TO JAEDA150 

IIG SHAPE/IIX.8F12.7) I AEDA160 

6A FORMAT l//10X,53HFLEXI 8LE SLOPE INCREMENTS AT LOAD POINTS DUE TO AAEOA170 
ILPHA/IIX,8F12.7) ) AEDA180 

65 FORMAT I//10X,53HFLEXieLE SLOPE INCREMENTS AT LOAD POINTS DUE TO 0AE0A190 

1ELTA/I1X,8F12.7) ) AEDA200 

66 FORMAT I//10X, 53HFLEXI BLE SLOPE INCREMENTS AT LOAD POINTS DUE TO AAEDA210 

IN, G/(1X,8F12.7) > AE0A220 

67 FORMAT I//5X,13HLOAD FACTOR =,F9.5) AE0A230 

68 FORMAT I /lOX ,1AHTRIM VALUES (,I3,12H I TERATIONS ) /3X, AHCA =,F10.6,AE0A2A0 
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APPENDIX — Continued 


13X,7HALPHA = ,F10 . 6 t 3X , 7HPH I = ,F10. 6/3 X,AHCN = tFlO.6 , 3X, 7HDELTA =AED4250 

2,F10.6,3X,7HAN, G =, F10.6/3X,4HCM =, F 10 .6 , 3X ,7HTHETA =, FIO . 6t 3X, 7HAED4260 
3QC/2V =fFlC.6) AE04270 

69 FORMAT ( /lOX , 24H0ER I V AT I VE CONTRIBUT IONS/8X, 9HJ 16 SHAPE 1 5X , 5HALPHAAED4280 

1,7X,5HDELTA,7X,5HQC/2V,7X,5HAN. G # 8X , 4HODQT , 8X t 4HMACH, 5X , 9H0rN PREAE04290 
2SS/3F CA,F26.7,5F12.7,F13.8/4H CN ,7F 12. 6,F 1 3. 7/4H CM , 7F1 2 . 6 , FI 3 . AE04300 
37) AE04310 

70 FORMAT < /lOX, 21HS TAB IL I TY OER I VATI VE S /12X , IHU ,10X ,4HUDOT, 7X , 5HALPHAE04 320 

1A,5X,9HALPHA DOT , 5X, SHTHETAf 9X , IHQ, 10X,4HODOT,7X , 5H0ELT A, 6Xt SHALT I AED433C 
2TUDE/5H CA , 8F1 2 .7, El 2. 5/4H CN 1 8F12 , 6 , E13. 5/4H CM , 8F12. 6 , E13. 5/ AE04340 
34H X ,8F12.6f E13.5/4H Z .8F12.6, E13.5/4H M i 8F12.6. E 13 . 5 ) AED4350 

71 FORMAT (/10X.17HSTATIC PARAMETERS/2C H CM ALPHA/CN ALPHA = , FIO .6/14AE04360 

IH STATIC MARGIN. 5X,1H=, FIO. 6/20H MANEUVER MARGIN =,F1C.6/8H DELTAED4370 
2A/L,11X,IH=,F10.6/12H CELTA/AN, G. 7X , 1H= , FIO . 6 ) AED4380 

72 FORMAT ( //lO X, 43HT0T AL TRIM CONDITION SLOPES AT SLOPE POINTS/I IX, 8AED4390 

1F12.7)) AED4400 

73 FORMAT ( //lOX, 42HT0TAL TRIM CONDITION SLOPES AT LOAD POINTS/ ( IX, 8FAED4410 

112.7)) AED4420 

74 FORMAT I //5X ,14HUNABLE TO TR I M/3X , 3HCN= , 8F12 . 6/3X , 3HCM= , 8F12 .6 » AED4430 

END AE04440- 
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APPENDIX - Continued 


I mil 


SUBROUTINE TRIM (CN*CM,W*FL#LFT,Q, $,Cf V,G, YI ,Y1 , Y2, AL 

,DEL, AN,CNT,CTRM 10 

lMTtQFrPR,PHIpTfIl) 

TRM 20 

OIHENSICN CN(8), CH(8I 

TRM 30 

WOOS=k/(Q*S) 

TRM 40 

G0G0=G/32.174 

TRM 50 

A=CM(2)/CM(3>-CN<2)/CN(3I 

TRM 60 

8= (WOQS-CN(5l l/CNI3l+CM(5)/CM(3l 

TRM 70 

IF (LFT.GE.l.ANO.FL.GT.GQGO) GO TO 3 

TRM 80 

QF=(32.174*FL-G)*C/(2.+V*V) 

TRM 90 

D=(CNJn+QF*CN{4) l/CNC 3 )•( CM ( 1 H-QF+CM ( 4 ) )/CM C 31 

TRM 100 

ALl=IO-FL*B)/A 

TRM 110 

DO 1 1=1,100 

TRM 120 

11=1 

TRM 130 

AN=FL*COS( ALl) 

TRM 140 

AL=(D-AN»B)/A 

TRM 150 

IF (AeS(AL-ALn*LE,*0000005l GO TO 2 

TRM 160 

AL1=AL 

TRM 170 

COATINUE 

TRM 180 

AN=FL*CCS< AL) 

TRM 190 

OEL=(-CM( n-AL*CM(2)-QF*CM«4)-AN*CM< 5)) /CM(3 ) 

TRM 200 

CNT=WOQS*AN 

TRM 210 

CMT = 0. 

TRM 220 

PR=0^ 

TRM 230 

PHI=0. 

TRM 240 

T=AL 

TRM 250 

RETURN 

TRM 260 

FW=WOCS*GOGO 

TRM 270 

FI=YI/(Q+S*C) 

TRM 280 

F1=FI*Y1 

TRM 290 

F2=FI*Y2 

TRM 300 

AL1=0. 

TRM 310 

TT=0. 

TRM 320 

CT=1. 

TRM 330 

sr=o. 

TRM 340 

CA = 1. 

TRM 350 

SA=0. 

TRM 360 

I=C 

TRM 370 

CP=GOGC/FL 

TRM 380 

QV=FL/GCGO-CT*CP*CA-ST*SA 

TRM 390 

Q1 = 0V*G/V. 

TRM 400 

QF=Q1*C/(2.»V) 

TRM 410 

QP=0V*CA/(2*»CT) 

TRM 420 

CP=SQRTa.-OV*TT*SA/Cr+GP*QP)-QP 

TRM 430 

IF (CP-1,) 6,5,5 

TRM 440 

R=C* 

TRM 450 

P = C. 

TRM 460 

SP = 0, 

TRM 470 

GO TO 7 

TRM 480 

SP=SQRT(l,-CP*CP) 

TRM 490 

P=-gi*TT/SP 

TRM 500 

R=G1*CP/SP 

TRM 510 

AN=GOGO*ICT*CP+QV*CA) 

TRM 520 

PR=P*R 

TRM 530 

CMT=F2*(P*P-R*R)-F1*PR 

TRM 540 

IF ( I.LT.O) GO TO 9 

TRM 550 

0=(CMT-CM( 1 )-QF*CM(4)+PR+CM( 6) )/CM(3 )+(CN( 1 )+QF*CN(4)- 

PR*CN(6) )/CNTRM 560 

1(3) 

TRM 570 

AL=(D-AN*B)/A 

TRM 580 

CA = COS(AU 

TRM 590 

SA=SIN( AL) 

TRM 600 
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APPENDIX — Concluded 


TT=SA*CP/CA 

TRM 

610 

T=ATAN2(SA*CPfCA> 

TRM 

620 

CT=COS(T) 

TRM 

630 

ST=SIN(T) 

TRM 

640 

IF ( ABSCAL-ALl) .LE.. 0000005) GO TO 8 

TRM 

650 

IF ( I.GE.lOO) GO TO 8 

TRM 

660 

1=1 + 1 

TRM 

670 

11=1 

TRM 

680 

AL1=AL 

TRM 

690 

GO TO 4 

TRM 

700 

I=-l 

TRM 

710 

GO TO 4 

TRM 

720 

CNT-WCQS*AN 

TRM 

730 

DEL=(CMT-CM( l)-AL*CM(2)-QF*CM(4)-AN+CM(5)+PR*CM(6) )/CM(3) 

TRM 

740 

PHI=ATAN2(SPfCP) 

TRM 

750 

RETURN 

TRM 

760 

END 

TRM 

77G 
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Figure 1. - Dlustration of reference axis system and several fundamental parameters. 
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Figure 2. - Panel arrangements used to represent example airplane. 






Dynamic pressure, q, kN/m^ Dynamic pressure, 5 1 kN/mz 


(b) M = 2.7. 

Figure 4.- Concluded. 
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(a) Concluded. 
Figure 5.- Continued. 






(b) M = 2.7. 
Figure 5.- Continued. 
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Figure 5.- Concluded. 
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(a) Concluded. 
Figure 6.- Continued. 
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(b) M = 2.7. 

Figvire 6.- Continued. 
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Figure 6.- Concluded. 
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(a) M = 0.8. 

Figure 7.- Partial derivatives of with respect to the physical variables. 
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(b) M = 2.7. 
Figure 7.- Continued. 



5( 




,m2/kN 




. 002 |- 


ko 



Weight, 2.67 MN 



10 20 30 40 50 

Dynamic pressure, g,kN/mz 


Weight, 1.87 MN 



0 10 20 30 4( 

Dynomic pressure ,q, kN/m^ 


(b) Concluded. 
Figure 7.- Concluded. 




Figure 8.- Variation of sonic velocity 
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(a) M = 0.8; c.g. at 0.41c. 

Figure 9. - Variation of characteristic roots with dynamic pressure. Rigid airplane. 
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(b) M = 0.8; c.g. at 0.47c. 
Figure 9.- Continued. 
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(c) M = 2.7; c.g. at 0.41c. 
Figure 9. - Continued. 
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Figure 9.- 
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(a) M = 0.8; c.g. at 0.41c. 

Figure 10. - Variation of characteristic roots with dynamic pressure. Elastic airplane. 
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(b) M = 0.8; c. 
Figure 10.- C 
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Figure 10.- Continued. 
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Figure 11.- Effect of root coupling on frequency and damping of short -period and phugoid modes. 
Elastic airplane; M = 0.8; c.g. at 0.41c; imiform atmosphere. 
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Figure 12.- Bode plots of frequency response amplitude of ujdp and np^6p. Elastic airplane 

M = 0.8; c.g. at 0.41c; weight, 2.67 MN. 
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(a) M = 0.8. 

Figure 13. - Variation of static control parameters with dynamic pressure. Rigid airplane. 
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(b) M = 2.7. 
Figure 14.- Concluded. 
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Figure 15.- Variation with dynamic pressure of approximations to static stability margins. 

Rigid airplane; uniform atmosphere. 
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(b) M = 2.7. 

Figure 16.- Concluded. 
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(a) M = 0.8. 

Figure 17.- Effect of dynamic -pressure contribution to stability derivatives with respect to u on 
characteristic roots. Elastic airplane; c.g. at 0.41c; uniform atmosphere. 
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Figure 17.- Concluded. 
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Figure 18.- Effect of dynamic -pressure contribution to stability derivatives with respect to u on static 
control parameters as calculated by using the present formulation and an indirect formulation. 
Elastic airplane; c.g. at 0.41c; uniform atmosphere. 
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(a) M = 0.8. 

Figure 19. - Comparison of the characteristic roots calculated by using the present formulation with those calculated 
by using an indirect formulation. Elastic airplane; c.g. at 0.41c; imiform atmosphere. 




